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Chapter 1 
Introduction 



1.1 Overview 

When we think of new physics beyond the standard model, supersymmetric (SUSY) 
extension |]T| of the standard model is one of the most attractive candidates. Cancellation 
of quadratic divergences in SUSY models naturally explains the stability of the electroweak 
scale against radiative corrections ^. Furthermore, if we assume the particle contents 
of the minimal SUSY standard model (MSSM), the three gauge coupling constants in the 
standard model meet at ~ lO^^GeV ^, which strongly supports grand unified theory 
(GUT) based on SUSY |, |^. 

In spite of these strong motivations, no direct evidence of SUSY (especially superpartners) 
has been discovered yet. Therefore, the SUSY is broken in nature, if it exists. Although 
many efforts have been made to understand the origin of the SUSY breaking, no convincing 
scenario of SUSY breaking has found yet. Nowadays, many people expect the existence of 
local SUSY {i.e. supergravity) and try to find a mechanism to break SUSY spontaneously 
in the framework. In the broken phase of the supergravity, super-Higgs effect occurs and 
gravitino, which is the superpartner of graviton, acquires mass by absorbing the Nambu- 
Goldstone fermion associated with SUSY breaking. In this case, the gravitino mass 1713/2 is 
expected to give us some informations about the SUSY breaking mechanism. For example, in 
models with the minimal kinetic term, the following (tree level) super-trace formula among 
the mass matrices Mj's holds; 

StrM2= ^ (-1)2^(2J+ l)trM^ - 2(r2^- l)m^/2, (1.1) 

spin J 

where is the number of the chiral multiplets in the spontaneously broken local SUSY 
model. In this case, all the SUSY breaking masses of squarks and sleptons are equal to 
the gravitino mass at the Planck scale. Meanwhile in models with "no-scale like" Kahler 
potential [1^ , SUSY breaking masses for sfermions vanish at the gravitational scale and 
are induced by radiative corrections, and hence the gravitino mass is not directly related to 
the scale of the SUSY breaking in the observable sector (which contains ordinary particles 
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in the standard model and their superpartners) . In order to understand the physics of the 
SUSY breaking, it is significant to clarify the property of the gravitino. But contrary to our 
theoretical interests, we have no hope to see the gravitino in collider experiments since its 
interaction is extremely weak. 

On the contrary, cosmological arguments provide us some informations about the grav- 
itino. In general, cosmology severely constrains properties of exotic particles. Let us review 
the constraints derived from cosmology. 

• The first is on the mass density of the exotic particle during the big-bang nucleosyn- 
thesis. If it is too large, it speeds up the expansion rate of the universe during that 
epoch and results in too many ^He. 

• The second is on the entropy production by the decay of the exotic particle. If the 
decay of the exotic particle releases a large amount of entropy, the baryon-to-photon 
ratio may become much well below what is observed today. 

• The third arises from the effects of the decay products on the big-bang nucleosynthesis. 
If the photon or some charged particle is produced by the decay of the exotic particle 
after the big-bang nucleosynthesis has started, energetic photons induced by the decay 
products may destruct light nuclei (D, ^He, '^He) and destroy the great success of the 
big-bang nucleosynthesis. 

• Furthermore, one can obtain the fourth constraints by considering the cosmic mi- 
crowave background distortion by the exotic particle with lifetime larger than ~ 
lO^^sec. 

• If exotic particle is stable, its present mass density provides us a fifth constraint. 

In fact, the most severe constraints on models based on supergravity are derived from the 
light element photo-dissociation and the present mass density of the universe. 

Following the above arguments, we can obtain stringent constraints on the gravitino mass 
in the standard big-bang cosmology. If the gravitino is unstable, it may decay after the big- 
bang nucleosynthesis and releases tremendous amount of entropy, which may confiict with 



the big-bang nucleosynthesis scenario. As Weinberg first pointed out [|Tl|], the gravitino mass 
should be larger than ~ lOTeV so that the gravitino can decay before the big-bang nucle- 
osynthesis starts. Furthermore in SUSY models with i?-parity invariance, unstable gravitino 
produces heavy stable particle {i.e. the lightest superparticle) in its decay processes, which 
results in unacceptably high mass density of the present universe |T^ . In order to reduce the 



number density of the lightest superparticle through pair annihilation processes, gravitinos 
should decay when the temperature of the universe is higher than (1 - 10)GeV. This requires 
that the gravitino mass should be larger than (10^ - 10'')GeV, which seems to be disfavored 
from the naturalness point of view. In the case of stable gravitino, the gravitino mass larger 
than ~ IkeV is excluded since the present mass density of the gravitino exceeds the critical 
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density of the present universe [O]. The above constraints on the gravitino mass seem to be 



very stringent especially for models with the minimal Kahler potential, since in such models 
the gravitino mass is expected to give the scale of the SUSY breaking in observable sector. 
In the inflationary universe W^, however, situation changes |TH|. In this case, the initial 



abundance of the gravitino is diluted during the inflation, and hence the number density of 
the gravitino becomes much less than that in the case of the standard big-bang cosmology. 
But even in the inflationary universe, the gravitino may cause the cosmological problems 
mentioned above since secondary gravitinos are produced through scattering processes off 
the background radiations or decay processes of superparticles. As we will see later, number 
density of the secondary gravitinos is approximately proportional to the reheating tempera- 
ture after the inflation and hence the upperbound on the reheating temperature is derived. 

In this thesis, we study details on the gravitino production in early universe and on its 
effects in the inflationary universe. Compared with the previous works, we have made an 
essential improvement on the following points. 

• Gravitino production cross sections are calculated by using full relevant terms in the 
local SUSY lagrangian. 

• High energy photon spectrum induced by the gravitino decay is obtained by solving 
the Boltzmann equations numerically. 

• Time evolutions of the light nuclei (D, ^He, ^He) with non-standard energetic photons 
are calculated by modifying Kawano's computer code. 

In our analysis, we assume that the light elements are synthesized through the (almost) 
standard scenario of the big-bang nucleosynthesis (with baryon-to-photon ratio 10~^ - 10"^''), 
and take the reheating temperature as a free parameter. 



1.2 Organization of this thesis 

The outline of this thesis is as follows. The former half of this thesis is devoted to 
the review of related topics, especially that of the gravitino properties. In Chapter ^ we 
review the motivation of SUSY. In Chapter ^, the gravitino field which is the gauge field 
associated with local SUSY transformation is introduced. Furthermore, lagrangian based 
on local SUSY is also shown in Chapter ^ and the super-trace formula in that framework is 
derived. Conventions used in Chapter ^ (and in other chapters) are shown in Appendix ^ 
In Chapter |^, we quantize a massive gravitino field and derive Feynman rules for gravitino. 

In the latter half of this thesis, we study the cosmology with the gravitino in detail. 
Overview of phenomenology with the gravitino is given in Chapter ^. In Chapter ^ ef- 
fects of unstable gravitino in the inflationary cosmology are analyzed in detail. In deriving 
constraints, we first derive photon spectrum induced by the decay of the gravitino. The 
procedure to obtain the photon spectrum is given in Appendix Then, we calculate the 
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time evolution of light nuclei with the obtained high energy photon spectrum, and we derive 
constraints on the reheating temperature and on the gravitino mass. In our analysis, we 
assume the standard big-bang nucleosynthesis scenario which is reviewed in Appendix |C[ 
The case of stable gravitino is discussed in Chapter [7[ Chapter |^ is devoted to discussions. 
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Chapter 2 

Motivations of supersymmetry 



2.1 Hierarchy problem in the standard model 

For particle physicists, symmetries in nature are significant guiding principles. Espe- 
cially, interactions of elementary particles (like quarks and leptons) can be understood by 
using the concept of the local gauge symmetry. Strong interaction is expected to originate to 
SU(3)(^ gauge group, and its theoretical predictions (like three gluon vertex and asymptotic 
free nature of its gauge coupling constant) have been confirmed experimentally. Meanwhile, 
results of recent electroweak precision measurements are in good agreements with the predic- 
tions of the spontaneously broken SU(2)^ x U(l)y gauge theory. Accompanied by theoretical 
and experimental successes, the standard model, based on the SU(3)^ x SU(2)^ x U(l)y 
gauge group, is regarded as the established one which describes particle interactions below 
the energy scale ~ lOOGeV. 

But once we look up high energy scale, one unpleasant problem, which is called hierarchy 
problem, appears in the standard model. In the standard model, existence of the elementary 
scalar boson, i.e. Higgs boson, is assumed in order to cause a spontaneous breaking of the 
gauge symmetry SU(2)^ x U(l)y — > U(l)g^. This is the origin of the hierarchy problem. As 
one can easily see, radiative corrections to the Higgs boson mass squared 5mjj are quadrat- 
ically divergent. Therefore, if one assumes the existence of the cut-off scale of the standard 
model AcuT at which the parameters in the standard model are set by a more fundamental 
theory. Sniff ~ O^aA^^j,) where a represents the coupling factor. The relation between the 
bare mass squared and the renormalized one m|^^ is written in the following way; 

m%^R = mH,B + Sm%. (2.1) 

In order to give the electroweak scale correctly, the renormalized mass squared m% ji should 
be O((100GeV)^). On the other hand, if we assume a larger value of Aqut, Smjj increases 
quadratically and a fine tuning of m%^g is needed so that the renormalized mass squared 
rrijj ji remains O((100GeV)^). For example, if we assume the cut-off scale of the standard 
model to be at the Planck scale ~ lO^^GeV, both m^ ff and Stjih are O((10^^GeV)^) for a 
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~ 0(1), and they should be chosen as 

+ ~ O((10i^GeV)2) - O((10^^GeV)2) ~ O((100GeV)'). (2.2) 

This is a terrific fine tuning. Therefore, if one assume that the cut-off scale of the standard 
model is much larger than the electroweak scale, we have to accept an unbelievable fine 
tuning of Higgs boson mass. This is the hierarchy problem. In fact, this problem stems from 
the fact that there is no symmetry which stabilizes the electroweak scale |T6|, [1^ . In order to 
solve this problem, we hope that some new physics (in other words, some new symmetry) in 
which quadratic divergences do not exist at all, appears at a energy scale O(100GeV — ITeV) 
and solve this difficulty. 



2.2 Supersymmetric extension of the standard model 

One of the most attractive solution to the hierarchy problem is SUSY SUSY is a 
symmetry which transforms bosons into fermions and vice versa. Therefore, in SUSY models 
the number of bosonic degrees of freedom is equal to that of fermionic ones. As we will see 
later, quadratic divergence of the Higgs (and other) boson masses are canceled out between 
the contributions from boson and fermion loops. 

Experimentally, however, we have not found any superpartners of the observed particles. 
This fact indicates that SUSY is broken in nature, if it exists. In order to solve the hierarchy 
problem, the SUSY must be broken softly |TB[ so that quadratic divergences do not exist at 
all. Usually, such a softly broken global SUSY model is regarded as a low energy effective 
theory of the spontaneously broken local SUSY model. We will comment on this point in the 
next chapter and here, we consider a phenomenologically acceptable (softly broken) SUSY 
model. 

When we extend the standard model to the supersymmetric one, we usually add "super- 
partners" for the ordinary particles existed in the standard model. In Table |2.1|, we show the 



particles in the minimal SUSY standard model (MSSM) and their gauge quantum numbers. 
Along with the existence of the superpartners, one big difference between the standard model 
and the SUSY one is the number of Higgs doublets, i.e. the MSSM requires two Higgs dou- 
blets (see Table p.l| ). In the SUSY standard model, Higgs bosons are accompanied by their 



fermionic superpartners which have the same gauge quantum numbers as the Higgs bosons. 
In this case, anomaly cancellation is not guaranteed if both of Hi and H2 are not included. 
Furthermore, in order to give fermion masses to up-type quarks as well as down-type quarks 
and leptons from Yukawa couplings of Higgs bosons, at least two chiral superfields Hi and 
H2 are needed. Mainly from the above two reasons, two Higgs doublets with representation 
(1, 2, -1/2) and (1, 2, 1/2) are introduced into the MSSM. 

Next, we will see the lagrangian of the MSSM. As a first step, we comment on im- 
parity. If we assume a particle content of the MSSM shown in Table p.l|, we can write down 
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Gauge sector 



Representation 


Boson {R = +1) Fermion {R = —1) 


(8,1,0) 
(1,3,0) 
(1,1,0) 


Gf, g 
w 

B, b 


Higgs sector 


Representation 


Boson {R = +1) Fermion {R = —1) 


(1,2,-1/2) 
(1,2,1/2) 


Hi xm 
H2 XH2 


Quark / lepton sector 


Representation 


Boson {R = —1) Fermion [R = +1) 


(3,2,1/6) 
(3*, 1,-2/3) 
(3*, 1,1/3) 
(1,2,-1/2) 
(1,2,1) 


Qi Qi 
ul ul 

# dl 



Table 2.1: Particle content of the minimal SUSY standard model. Index i is the generation 
index which runs from 1 to 3. For each particles, representation of the SU(3)(^ x SU(2)^ x 
U(l)y gauge group is also shown. 



interactions which violate baryon- or lepton-number conservations. For example, interactions 
such as u'^df^df^ or d'^ql cannot be forbidden by gauge invariance or renormalizability. But 
phenomenologically, strength of these interactions is severely constrained since they may 
induce unwantedly high rate of nucleon decay and neutron-anti-neutron oscillation, and 
they wash out baryon number in the early universe |T3. Rather than assuming extremely 



small coupling constants for them, we usually forbid these dangerous terms by introducing 
a discrete symmetry, that is called i?-parity. i?-parity assigns -|-1 for ordinary particles 
in the standard model and —1 for their superpartners. One can see that if we require 
the invariance under the i?-parity, baryon- and lepton-numbers are conserved under the 
restriction of renormalizability. In this thesis, we adopt the i?-invariance below. Notice that 
i?-invariance also guarantees the stability of the lightest R-odd particle, i.e. the lightest 
superparticle (LSP). 

Assuming the i?- invariance, the superpotential of the MSSM is given by 

Wmssm = ylj^u'rqjH2 + ylfdlqjHi + y\felljHi 

+IIhHiH2, (2.3) 

where i and j are generation indices, and we have omitted the group indices for simplicity. 
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Here, ?/*^"\ y^'^^ and y^^^ are the Yukawa couphng constants of the up-, down- and lepton- 
sector, respectively. 

Since the SUSY should be broken in nature, SUSY breaking terms are also necessary in 
lagrangian. In order not to induce quadratic divergences, SUSY should be broken softly. In 
general, soft SUSY breaking terms are gaugino mass terms, scalar mass terms and trilinear 
coupling terms for scalar bosons of same chirality In the MSSM, SUSY breaking terms 
are given by 



C 



soft 



dij 



-mjj^ l-f^il^ - "^^2 1^2^ - (rnlHiH2 + h.c. 
- {jTLG^gg + mG2WW + mcihl + h.c. 



[2.A) 



where m? - m| are the squark and slepton masses, and mcz ^ the gauge fermion masses. 
In the next chapter, we will see that these SUSY breaking parameters can be obtained in a 
low energy effective theory of local SUSY models. 

As mentioned before, quadratic divergence of the two point functions of scalar bosons 
disappears in SUSY models. At the one loop level, this can be easily seen. For example, 
we will see the cancellation in the mass of H2. Feynman diagrams which give quadratically 



divergent radiative corrections to the H2 mass are shown in Fig. |2.1j . Each of them are 
quadratically divergent; 



CB 



CF 



GB 



GF 



in 



2 J Vij 



CUT 



+ 



47r2 V2 
An' 



^gi + -gt]AcuT + 



(2.5) 
(2.6) 
(2.7) 
(2.8) 



where CF (CB, GB, GF) represents the contribution from chiral fermion (chiral boson, gauge 
boson, gauge fermion), ■■■ the terms which do not contain quadratic divergences and Aqut the 
cut-off. These quadratic divergences, however, cancel out between boson and fermion loops. 
Quadratic divergences in other scalar masses also disappear in the same way, and hence the 
hierarchy problem can be solved by extending the standard model to the supersymmetric 



one. 



In the MSSM, other interesting new physics, i.e. the grand unified theory (GUT) pO| , 
is suggested from the renormalization group analysis 0, As mentioned before, SUSY 
extension of the standard model increases the number of particles, which changes the renor- 
malization group equation of the gauge coupling constants. In Fig. |2.2|, we show the renor- 
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Figure 2.1: Quadratically divergent Feynman diagrams for the if2 mass. Diagram with CF 
(CB, GB, GF) is the contribution from chiral fermion (chiral boson, gauge boson, gauge 
fermion) loop. Dashed hues in external lines represent H^- Notice that the diagram with H2 
loop (lower-left) originates to the gauge D-term, and hence we classify it as the contribution 
of gauge boson. 

malization group flow of SU(3)^, SU(2)^ and U(l)y gauge coupling constants in the MSSM 
case and in the standard model case. In the MSSM case, three gauge coupling constants 
meet at the energy scale ~ lO^^GeV which may be identified with the GUT scale, while in 
the standard model case, the renormalization group flow of the gauge coupling constants 
conflicts with the gauge coupling unification. 

In 



Another indirect evidence of SUSY GUT is the bottom-tau mass ratio 21, 22, 23 



SU(5) or SO(IO) GUT, Yukawa coupling constants of down- and lepton-sectors to the Higgs 
boson are also expected to be unified at the GUT scale, and hence we can get a rela- 
tion between the bottom- and tau- Yukawa coupling constants at the electroweak scale. By 
using this relation, the mass of the bottom quark is obtained once the tau-lepton mass 
(^ 1.777GeV |4 
function of tan /3 



is fixed. In Fig. |2.3| , we show the predicted bottom-quark mass as a 
(iJi) / (-^2)- Notice that the maximal and minimal values for tan/? are 
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10' 10^ 10^ 10^ 10^ 10^ 10' 10^ 10^10'°10''10'^10'^10'^10'^10'^10''10'^ 

(GeV) 

Figure 2.2: Renormalization group flow of tlie coupling constants of SU(3)p, S\J{2)j^ and 
U(l)y gauge group for the case of (a) the MSSM, and (b) the standard model. Here, we use 
two loop renormalization group equations, and take the SUSY scale at ITeV for the MSSM 
case. 
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Figure 2.3: The predicted value of the running bottom-quark mass mf,(mi,) is shown as a 
function of tan/3 for as^niz) = 0.11 and 0.12. Here, we take the (on-shell) top-quark mass 
at 174GeV. 



determined so that all the Yukawa coupling constants do not blow up below the GUT scale. 
As one can see, SUSY GUT predicts the bottom-quark mass to be (4 - 6)GeV which is 
close to that determined from experiments; mf,(mf,) = (4.25 ± 0.2)GeV |^ (where we have 
doubled the uncertainty), especially when tan/5 approaches its maximal or minimal value. 

Contrary to those attractive features, no direct evidence of SUSY has been found, which 
certainly indicates that the SUSY is a (softly) broken symmetry. The physics of SUSY 
breaking is, however, still an open question and we have not understood it yet. Especially in 
the framework of global SUSY, it seems to be very much difficult to construct a phenomeno- 
logically favorable model. One of the reason is that there exists a mass formula in the global 
SUSY model; 



StrM' 



j=o 



'trM^ 



0, 



(2.9) 



which prevents all the squarks and sleptons from having masses larger than those of quarks 
and leptons. To avoid this constraint, many people extend the global SUSY to the local one 
and consider the physics of SUSY breaking in the framework of supergravity. In the next 
chapter, we will investigate local SUSY model and see how the mass formula in global SUSY 



models ( |2.9| ) is modified. 
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Chapter 3 

Review of supergravity 



In this chapter, we will introduce a lagrangian which is invariant under the local SUSY 
transformation, and derive super-trace formula in that framework. Conventions used in this 
chapter are essentially equal to those used in ref.[^ except that we use the metric (in flat 
space-time) as g^j^y ~ diag(l, —1, —1, —1). For our convention, see also Appendix ^ 



3.1 Heuristic approach to supergravity lagrangian 

Compared with the global SUSY, one of the characteristics of the local one is the 
existence of a gauge field associated with the local SUSY, which is called gravitino. As in 
the case of ordinary gauge theories, the gravitino couples to a Noether current of SUSY 
and maintains the invariance under the local SUSY transformation. In this section, we will 
briefly review the role of the gravitino in the local SUSY theory by using the simplest model, 
which is the Wess-Zumino model without interactions. 

Let us begin with the global case. In this case, total lagrangian contains only two terms, 
one is the kinetic term of a massless complex scalar boson and the other is that of a 
massless chiral fermion x'l 

Cy,^ = dyd^ct)* + ixcy>'d,x- (3.1) 

Up to total derivative, this lagrangian is invariant under the following global SUSY trans- 
formation, 

6<P = V2^x, (3.2) 
6x = -zV2a^ad,^), (3.3) 

where ^ is the infinitesimal Grassmann-odd parameter. 



If ^ has space-time dependence, lagrangian (|3.1| ) is not invariant but extra terms which 
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are proportional to dC, or dC, appear with the supertransformation (|3.2|) and 



5£wz = V2{{d^0(^''cx^'x{du<P*) + w''(r''{d^O{du<P)} + (total derivative) 

= i (9^0 + ^-c- + (total derivative). (3.4) 

where 

Jf" = -iV2a''a''x{d^<l)*)- (3.5) 

Notice that is the Noether current of SUSY, which is called supercurrent. In order to 
keep invariance, we introduce a gauge field '^p^i■ As in the cases of ordinary gauge theories, 
the gauge field couples to the supercurrent in the following way; 

C^j = -'-GsiJ^J^ + h.c. , (3.6) 

where Gs is the "coupling constant" which we will determine later. Since the charge of SUSY 
has a Grassmann-odd nature with spin index, the gauge field ^/'^ associated with SUSY is a 
spin I fermion. Varying eq.( p.6|) , one obtains 

SC^j = -'-Gs {iSij^)J^ + i^f^iSJ^} + h.c. (3.7) 



Therefore, if ipu, transforms as 



2 

Sipt^ ~ -^d^^, (3. 



the first term in eq.(|3.7|) cancels out the contribution from eq. 



Next, we will consider the second term in eq. ( ^.7] ). Supertransformation of the supercur- 
rent ( |3.5| ) gives energy-momentum tensor T^" of the chiral multiplet {4>,x)'i 

{Qa^Ja} = -2<^T/ + (total derivative), 

{Q",T"} = -2a'^'^"T/ + (total derivative), (3.9) 

where Q and Q are the generators of the SUSY transformation, and hence the second term 
in eq. (|3.7|) becomes 

- ^Gs {M^Jn} + h.c. = '-Gs {^.aj^ + ^.a^? + + T^^ (3.10) 

In order to cancel out these terms, we rewrite the lagrangian ( |3.1|) by explicitly expressing 
the metric tensor g^^; 

jOwz ^ V^g,. {d^(Pd-(P* + txcy^d'x) , (3.ii) 
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where g = det g^,y, and use the fact that the metric is the gauge field associated with the 
energy- momentum tensor (which is the Noether current of the space-time translation), that 
is, the energy-momentum tensor T^^ is obtained if one varies lagrangian by (yf^j^; 

dC 

- 2T^^ (3.12) 



dg, 



Then, the metric tensor g^^, {i.e. graviton) can be regarded as a superpartner of the gravitino 
field ipiJ., and its transformation law is determined so that the local SUSY invariance is 
maintained; 

~ -^Gs {rcr"^ + r^^l + ^'o'^i + W'^'i) T,u. (3.13) 



Combining eq.( |3.12| ) with eg. ( ^.131 ), one can obtain the transformation law of the metric 
tensor; 

^g^lv iGs (^^(x^? + i^vfTi^i + i^^^^ui + ■ (3.14) 

In the following arguments, in fact, it is more convenient to use the vierbein e^" rather 
than the metric tensor g^^, = r]abe-y°'ej' , where r]ah = diag(l, —1, —1, —1) is the metric tensor 
in flat space-time. In supergravity models, transformation law of the vierbein is deflned 
as 

<5e/ = -zGs {^a^i^, + e^>^) , (3.15) 



with Gs = M ^. This transformation law gives eq.( |3.14D . Under the local SUSY transfor- 



mation, the vierbein e^" and the gravitino ip^ (and some other auxiliary fields) make up a 
multiplet, which we call a supergravity multiplet. 

As we have seen, if we extend the global SUSY to the local one, the metric tensor g^j^^, 
automatically comes into the theory, and hence we must consider gravity. This is the reason 
why the local SUSY is sometimes called supergravity. 



3.2 Minimal supergravity model 

In the previous section, we have introduced the gravitino field in order to keep the 
local SUSY invariance. As we have seen, the gravitino field couples to the supercurrent, but 
the strength Gs of the coupling between the gravitino and the supercurrent has not been 
determined yet. In order to determine the coupling strength G^, we must see the invariance 
of the kinetic terms of and under the local SUSY transformation. 

In this section, we will explicitly investigate the local SUSY invariance of the minimal 
supergravity model which contains only the graviton and the gravitino field ipfi- As a 
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result, we will see the coupling strength Gs should be equal to the inverse of the gravitational 
scale. 

The lagrangian of the minimal supergravity model is given by 



-^MSG — -^EH + ^RS, (3.16) 

where £eh is the usual Einstein-Hilbert lagrangian, and £rs is the Rarita-Schwinger la- 
grangian which is essentially the kinetic term of the spin | gravitino field. The Einstein- 
Hilbert lagrangian is given by 

Ceh = -— ei?, (3.17) 

with 

R = ea^et" (9^0;/' - d^u^'^' - u^^'uj + ^^^^c) , (3.18) 

where cu^"^ denotes the spin connection and e = det Cq'^. On the other hand, the Rarita- 
Schwinger lagrangian can be written as 

^RS = ee^^P^'^^a.'bp^,, (3.19) 

where 6^'^/'°" is the totally anti-symmetric tensor (eoi23 = — 1 in flat space-time) and the 
covariant derivative of the gravitino field is given by 

V^i,, = d^ij, + Kj^^^aabi'u. (3.20) 



In the following, we will see the invariance of the minimal supergravity lagrangian ( p. 161) 
under the local SUSY transformation; 

= -iGs + la'^iJ,) , (3.21) 



= = 1^ {d,i + ^cu/V,,e) , (3.22) 

where parameter Gs will be determined so that the local SUSY invariance is maintained 
(see eg. ( |3.8| ) and eg. ( |3.15|) ). 



Before checking the invariance of the minimal supergravity lagrangian ( p.l6| ), we will 
comment on the spin connection uj^""^ . As we will see below, uj^'^^ is represented as a function 
of the vierbein e^" and the gravitino ip^ by solving its field eguation; 



15 



and hence uo^""^ is regarded as an auxiliary field. From eq.( |3.23D , one can obtain 



(3.24) 



By solving this equation, explicit form of the spin connection is given by 



ab 



77 (^M^P dp^fj, ) ~l~ ^pa iS^u^p d^Cfj ) a pa idpC^j f^cr^p ) }" 



(3.25) 



Now let us see the invariance of the lagrangian ( p.l6| ). With the help of chain rule, 
variation of the total lagrangian is given by 



5C 



MSG 



5e 



MSG 



MSG 



(3.26) 



where we have used the fact that the spin connection u;^"^ is a function of the vierbein e^"' 
and the gravitino ipp- The important point is that the last two terms in eq. (|3.26|) which 
are proportional to (5£msg/^^^) vanish since the spin connection obeys its field equation 
{SCmsg/S^^) = 0. Therefore, we only have to vary e^"" and ipp (with u^'^^ fixed) in order to 
obtain 5£msg- (We denote this operation A.) 

Variation of the Einstein-Hilbert lagrangian £eh gives the Einstein tensor multiplied by 
"^P ) 



AC 



EH 



6C 



EH 



iGsM'e [ - -ea^R]ij^a-^ + h.c. 



(3.27) 



with 



Rp"^ = ^b 



ab 



ab 



ac, , b 



dv^p"" ~ ^p""^^vc 



^v°''^^pc 



(3.28) 



On the other hand, after some straightforward calculations, A£rs becomes the following 
form; 



RS 



r a 



5C 



RS 
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-ee^"""" {J- (5^e/ + u;/^e,f.) + (V^^cr>,) | laaVp^,. (3.29) 



By setting Gs = M ^, the first line in eg. ( p. 291) is equal to —ACer, the second line vanishes 



due to eq.( p.24|) , and hence 5£msg vanishes; 

SCmsg = A£eh + AC^s = 0. (3.30) 

This is the end of the proof of the invariance. We have seen the invariance of the minimal 
supergravity lagrangian ( p. 161) under the local SUSY transformation (|3.21|) and (|3.22|) with 

Gs = M-\ 

This fact suggests that the coupling strength Gs between the gravitino field and 
the supercurrent is not a free parameter but a model independent constant which is 
determined by the requirement that the kinetic term of the supergravity multiplet {ip^, e^") 
is invariant under the local SUSY transformation. In the next section, we can see that 
general supergravity lagrangian contains interaction terms between the gravitino field ipfj, 
and the supercurrent with definite couphng strength {i.e. Gs = M~^)] 

^^J = -^^Pf^J' + ^-c- (3-31) 

As we will see in the following chapters, such interaction terms become very important in 
investigating phenomenology with the gravitino. 

3.3 General supergravity lagrangian 

In this section, we will extend the minimal supergravity model to general one. Derivation 
of the general supergravity lagrangian is given in elsewhere [B], |2^, but it is very much 



complicated task. Therefore, we only give a final form here by following ref.|^. In this 
section, we use the M = 1 unit for simplicity. 

The general supergravity lagrangian is essentially characterized by three functions; Kahler 
potential K{(l),(f)*), superpotential W{(f)), and kinetic function /(0) for vector multiplets. 
Notice that the Kahler potential K{(j),(l)*) is a function of scalar fields and 0*, while the 
superpotential 11^(0) and the kinetic function /(0) depend scalar fields with definite chirality. 

By using these functions, the general form of the supergravity lagrangian, which contains 
scalar field 0, chiral fermion x, gauge boson A^, and gauge fermion A as well as the vierbein 
Cfj^"" and the gravitino ip^ can be written as 

/:suGRA = -^eR + egi,,V,<p'V^<p*^ -^eg^D^,^D<^^^ 
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I 

+2^ 



7(a) 



(a)C 



(6) 



+V2e5ft,-*XJ)X^A(«) + V2e5^,,-*X('„)X'A 
-^V2ea,/(„,)xV'^^A(»)FW - iV2ea,,/(;,^x^a'^'^A^"VW 



2 

i 



-^e [c/ij.^w* - 2Rij*ki'] x'xVx^ 
1 

+ 16'= 



_|_ pi") 
up ' up 



XCTaX 



^9irftai>) + f''^'"'^''9ifii,c)dj.f^ x^a'^x^A^^^^.A^^) 
+^eV,a,/(„6)X^X^A('^)A('') + ^eV^.a^./f^.^^A^-^A^"^ 
+^e/«(^'^)-i9,/(„e)a,/(,,)X^A(«)x^A(^) 



-^e^^^*a./(,5)a,./(;,)A('^)A(^)A^^^A^'^ 

+ Y^eA(a)C7''A^"^A(6)(7^A^''^ 

+ '-V2edd\,,) [xV'^'^A^^Vm^.A^" - ^V^;.a'^x'A(»)A(^) 



c/^^* (AW") {Dj*W*) - 3W*W 



(3.32) 



where = Re/ and = Im/. Indices i, j, • • • represent species of chiral multiplets, and 
(a), (&),•■■ ^re indices for adjoint representation of gauge group (with gauge couphng constant 
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g and structure constant / ) which are raised and lowered with /^^^^^ and its inverse. Notice 
that the Kahler potential K and the superpotential W in the total lagrangian ( p.32|) can be 
arranged into the following form; 



G = K + \Yi{W*W) 



Covariant derivatives are defined as 



V 



F'^ is the field strength tensor for the gauge boson A'^^\ and is defined as 



Differentiation by the scalar field is symbolically represented by the index i; 



OA- ■ 



OA---) 



(3.33) 



(3.34) 



(3.35) 



+ - {kP.cI^^ - K,.V,<f>*^) A^-^) + ^^74'')lmF(,)A('^), (3.36) 



(3.37) 



(3.38) 



(3.39) 



With the help of eq.( 3.39| ), "derivatives" of the superpotential are defined as 



DiW = Wi + KiW, 



(3.40) 



ViDjW = Wij + KijW + KiDjW + KjDiW-KiKjW-T^jDkW. (3.41) 

"Metric of the Kahler manifold" gij* is defined by varying the Kahler potential K by the 
scalar fields 0* and (j)*^; 



gij* 



(3.42) 
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and g^^* is its inverse, 

g,,.g''* =6^:, g,,.g'^* = 6l. (3.43) 
From this metric, "connection" F^^ and "curvature" Rij*ki* is given by 

r& - (3-44) 
- ^^^3'^* - 3"""* {-ir^-^^) {^^^^^ (3.45) 
By using the connection given in eq. (|3.44| ), "covariant derivative" Vj is defined as 

V.VS- ^ - Ff^.\4. (3.46) 

(Here, is a function of scalar fields with index i.) 

Next, we will comment on X^""^^ and D^'^\ X*^°)*(0) is the Killing vector associated with 
Kahler metric gij*. That is, with the field transformation 

0* ^ 0*' = 0* + X(")^(0)e, (3.47) 
X*(")*((/)*)e, (3.48) 



(where e is a infinitesimal parameter), the Lie derivative 6^x^ of g^j* vanishes; 



= 0, (3.49) 

where xi"'' = gijX^"'^^ , and the index i represents both i and i*. From the above equation, 
we obtain two equations for the Killing vectors X*^"^*(0) and X**^")*(0*); 

V,xf ^ + VjXf"^ = 0, (3.50) 
ViXj?^ + Vj.xf"^ = 0. (3.51) 

The former equation is automatically satisfied, while the latter allows one to write down the 
Killing vectors X'^"^*(0) and X*(")*((/)*) as derivatives of some function D^"''^ (which is called 
Killing potential); 

XW^(0) = (3.52) 
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(3.53) 



By solving eq.(^3TD, eq. (|332|) and eq.(^33|), the Killing vectors X^'''^'{4>), X*^"'^' {(p*) , and the 
Killing potential D^"-* can be obtained. F^"'\ which is a analytic function of (f)\ is defined as 



For example for the minimal Kahler potential K 
Killing potential take the following forms; 



(3.54) 

the Killing vectors and the 

(3.55) 
(3.56) 
(3.57) 



where T^" is a generator of gauged Lie group. 

For the local SUSY transformation, variations of each component field are given by 



(3.58) 
(3.59) 



X 



2V2 



2V2' 



with 



(3.60) 
(3.61) 

(3.62) 
(3.63) 

(3.64) 



Notice that given in eq. (|3.64| ) corresponds to the auxiliary field in chiral multiplet in the 
global SUSY case. (Do not confuse with F^^^ defined in eq.( 3.54 ).) 
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As mentioned in the previous section, the total lagrangian ( ^.321 ) contains interaction 
terms between the gravitino field and the supercurrent in the form of eq.( p.31D . For a 
later convenience, we note the interaction terms here; 



e {^.cr'^Pa^a) + ^.a'^'^a'^Xia)) FI% (3.65) 



2M 

where we have explicitly written down the M-dependence. 

3.4 Super-Higgs mechanism 



With the general supergravity lagrangian ( ^.321 ), we now can investigate the super- 



Higgs mechanism |§] which is closely related to the spontaneous breaking of the local SUSY. 
If the global SUSY is broken spontaneously, massless Nambu-Goldstone fermion which is 
called goldstino appears, while in spontaneously broken local SUSY models, this goldstino 
component is absorbed by the gravitino through the super-Higgs mechanism. In this section, 
we will discuss this mechanism in more detail. (As in the previous section, we take M = 1 
unit in this and the next sections.) 

We will begin by considering the SUSY breaking in supergravity. In global SUSY mod- 
els, a spontaneous breaking occurs if some auxiliary field, which is obtained by the SUSY 
transformation of some chiral fermion x or gauge fermion A, receives non- vanishing vacuum- 
expectation values. In the local SUSY case, we also use 6x and 6X as order parameters 
of the SUSY breaking from the reasons mentioned below. As one can see in eq.( |3.6(J| ) and 



eq. p.62| ), the local SUSY transformations of spin | fermions contain terms which are similar 



to the auxiliary fields obtained in the global case. Assuming (local) Lorentz invariance of 
the ground state and no fermion- fermion condensation (like gaugino-gaugino condensation), 
{6x) 7^ and (SX) ^ reduce to the following conditions; 

{bx') = -V2{F^)i + 0, (3.66) 
(^A^'^)) = -z^(Z}("))^ ^ 0. (3.67) 

As we will see later, if 5% or bX has non-vanishing vacuum-expectation value, the gravitino, 
which is the gauge field associated with the local SUSY, absorbs a massless eigenstate of 
fermion mass matrix and acquires a non- vanishing mass (which is called super-Higgs mecha- 
nism). Furthermore, the vacuum-expectation value of or D'^"-^ may provide a mass splitting 
of bosonic and fermionic states [i.e. StrM^ ^ 0, see the next section). These phenomena 
indicate a spontaneous breaking of SUSY and hence in this thesis, we regard bx and bX as 
order parameters of SUSY even in local SUSY case. 

Before investigating detail of the fermion mass matrix, we give some comments. The 
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first comment is on kinetic terms of cliiral and vector multiplets. In supergravity lagrangian 
( |3.32|) , tliese kinetic terms are cliaracterized by tlie Kahler potential K and tlie kinetic 



function /, respectively. The kinetic terms are properly normalized if gtj* = Kij* = Sij* H , 

and f(ab) = Sab + ■ ■ ■, where ■ ■ ■ represents the higher order terms which produce interactions 
of higher dimensions. For simplicity, we ignore these higher order terms and use the minimal 
Kahler potential iiT™" and the minimal kinetic function J™"; 

ir-- = ^0>«, fl2l^ = 5ab. (3.68) 

i 

Another comment is on the cosmological constant. Since we are interested in field theories 
in the Minkowski space-time, we require the cosmological constant to vanish. In the super- 
gravity models, the vanishing cosmological constant is obtained if the following condition is 
satisfied; 

{{D,W) {Di,W*) - 3W*W} + ^^g"^ D^"^ D^"^ = 0. (3.69) 

If e^/^DiW or has a non-vanishing vacuum-expectation value, the superpotential W 
should have a non-vanishing vacuum-expectation value in order to satisfy the condition 
( p.69| ). As we will see later, the gravitino mass is given by e^^'^W, and hence the condition 



for the vanishing cosmological constant ( p.69 ) requires non- vanishing gravitino mass. 



Now, let us discuss the super-Higgs mechanism. We first write down the quadratic part of 
the fermion terms without derivative in the lagrangian (|3.32|) by using G defined in eq. ( |3.33|) ; 



+V2g (-iDfx'X^"^ + iDPt^^""^] . (3.70) 



Notice that if e^/^DiW or D^"^ has a non-vanishing vacuum-expectation value, the gravitino 
field ip^ mixes with spin ^ fermion x or A. These mixing terms can be eliminated by a shift 
of the gravitino field ipii] 



6 ^ 
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+t {y/2gDP - ^gG^.D^^^^ x^A^^^ + h.c. , (3.71) 

with 

V^^G,,t + \e-^''gD^''^\^''\ (3.72) 

From the above lagrangian, we can read off the gravitino mass ^3/2 asQ 

rri3/2 = e^/2^ |e^/2^|. (3.73) 

On the other hand, masses of the spin | fermions (x*, A^"-*) can be obtained from the following 
mass matrix 



(3.74) 



(1/2) (1/2) I , 
\ "^ja ^ab J 

whose each components are given by 

mSf) = (g,, + i^.G,) , (3.75) 

mlf) ^ le-^/2/D(")D(^), (3.76) 
miy^) = i^2[-gDf^^\gG,D^^^\^. (3.77) 



One can easily check that the fermionic field rj defined in eq.( |3.72|) is a massless eigenstate 



of the mass matrix ( 3.74| ) and hence it is a goldstino component. Therefore, the gravitino 



field -ip^ absorbs the goldstino component r] and acquires a non- vanishing mass m3/2 = e^^'^. 

3.5 Super-trace formula in supergravity 

Having seen the structure of the mass matrix of fermionic sector in the previous section, 
we will next investigate mass differences between bosonic and fermionic states. For this 
purpose, we use the super-trace of the mass-squared matrices 



3/2 

J=0 



StrM^ = ^(-l)^trM^, (3.78) 



since it represents some information on the mass splitting of bosons and fermions. 

First, let us consider trace of the scalar mass matrix. In deriving masses of scalar bosons 
(and other fields), we expand the scalar potential at a stationary point in scalar field. 

^For simplicity, in the following arguments in this and the next sections, we omit the bracket (O) which 
represents the vacuum-expectation value of O. 
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From the supergravity lagrangian (|3.32|), scalar potential V is given by 



V{<P, 0*) = {G,G,, - 3) + ^g^D^^^D^^\ (3.79) 
Stationary condition can be written as 

n - — 

= Gje^ {GiGi. - 3) + (GijGi, + GiGi.,) + g^D^'^^Df. (3.80) 

Multiplying this equation by Gj., and using the identity of^Gj* = one can obtain 

e^GjGj, {GiG^* - 3) + {GijGr*Gj, + GiG^*) + g^D^^^D^"^ = 0. (3.81) 

Furthermore, we demand the cosmo logical constant to vanish at the stationary point. The 
condition for vanishing cosmological constant has been derived in eq.( |3.69| ), which becomes 

{GiG^, - 3) + = 0. (3.82) 

The mass-squared matrix for scalar fields can be obtained from the second derivative of 
the scalar potential V. Especially, its diagonal part is given by 



dcjy^dcf)*^ — e'^|GjG'fc* {GiGi* — 3) + 6jk* {GiGi* — 3) + GijGi 
+Gj (Gfe. + GiGi^k") + Gk* (GijGi* + Gj) + 5jk*^ 
+g\Df^D^} + D^''^Df^,). (3.83) 



With the stationary condition ( |3.81D and the condition for vanishing cosmological constant 



( |3.82D , the trace of the scalar mass matrix Mq can be obtained as 



2 



+2e^ {G,,G,,,, + n^) + 2g' (ot^ + D^") A?) , (3.84) 



where = X)? Qa* is the number of chiral multiplets. 

Next, we will consider the mass matrix for vector bosons. Mass terms of vector bosons 
come from the covariant derivatives of scalar fields; 

^m5"''^4"M(^)'^ = ^V'7;"Tj'fc0''A(fU(^)^. (3.85) 
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From the above equation, one can easily read off the mass-squared matrix for vector bosons, 
and its trace is given by 



trM^ = 3 X 2/0*'Ti«r;^0'= 



(3.86) 



where the Kilhng potential D^'^^ for the case of the minimal Kahler potential is given in 
eq.( |3.57| ). Notice that the prefactor 3 in the right-hand side of eq.( |3.86|) corresponds to the 
number of polarization of a massive vector boson. 

The mass matrix of fermionic sector has been derived in the previous section, and we can 
easily obtain the trace of it. Masses of spin | fermions are obtained from the matrix ( |3.74| ), 
and the trace of the spin ^ fermion mass-squared matrix is given by 



trM^ 



/2 



2x E 



(1/2) 



+ E 



(1/2) 
ab 



ab 



2E 



(1/2) 



2e^ {G.jGi.j* - 1) - 2g^D^^'^D^'''^ 



-ie-^ ((^^z^C-) + 8g'Dl^'>Dl^\ 



(3.87) 



In supergravity, spin | particle is only the gravitino, and the trace of its mass matrix is 



given by 



trM 



3/2 



4e 



G 



Am 



3/2- 



(3. 



Combining eq. p.84| ), eq. ( p.86D , eq.( p.87D and eq.( |3.88D , the super-trace of the mass- 
squared matrices are given by 



StrM' = 2{n^- l)ml^ - {n^ - 1) g"" D^""^ D^""^ + 2g^D^^^DlP. 



(3.89) 



In SUSY models with e^^'^DiW ^ and D^"^ = {i.e. if the SUSY is broken by the F-term 
condensation), the SUSY breaking is characterized by gravitino mass m3/2, and scalar par- 
ticles are expected to become heavier than their superparticles. This is phenomenologically 
favorable. In the next section, we will see an example for such models. 



3.6 Model with Polonyi's superpotential 

In constructing a phenomenologically acceptable model based on supergravity, we mostly 
consider a model which contains two sectors; a so-called hidden sector responsible for the 
spontaneous breaking of SUSY, and an observable sector which contains ordinary particles 
such as quarks, leptons, Higgses, gauge fields, and their superpartners. The strength of inter- 
actions between the hidden and the observable sectors are at the same order of gravitational 
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one, and hence these two sectors couple very weakly for each other. In this section, we will 
investigate a simple model for the hidden sector proposed by Polonyi As we will see 
below, this model is simple but very suggestive. 

We first discuss the hidden sector. The simplest hidden sector, proposed by Polonyi, 
contains only one chiral multiplet (0p,Xp)) which takes the following superpotential Wp; 



(3.90) 



where /i and uj are the free parameters which we will determine by the phenomenological 
requirements. For the Kahler potential for the Polonyi field 0p, we use the minimum one, 
i.e. K = (j)*p4>p. In this model, the order parameter of SUSY breaking is given by 



M2 



+ iu) + l}e 



(3.91) 



For the case |ti;| < 2M, solution to the equation Fp = does not exist, and the SUSY 
is expected to be spontaneously broken. The scalar potential for the Polonyi field 0p is 
obtained as 



1 



bp + u) + l 



M2 



^P + cone'^-^p/^''. 



(3.92) 



The minimum of this potential is given by the solution to the following equation; 
dV 







M2 



{m2' 



Jp + Uj)\—(pp {(f)p + Uj) + l 



+<Pp[^<Pp {<Pp + Uj) + l 



+ UJ) 



M2 



M2 



+ UJ) + 1 



M2 



')p + UJ 



(3.93) 



Furthermore, since we are interested in field theories in fiat space-time, we demand the 
cosmological constant to vanish at the potential minimum; 



(V) 



M2 



+ UJ) 



M2 



0. 



(3.94) 



The parameter u is chosen so that eq.( ^.93| ) and eq. ( |3.94| ) have a solution simultaneously. 
Then cu is determined to be = ±(2 — \/3)M.0 Hereafter, we use the branch u = {2-V3)M. 



^In fact, even for uj = ±(2 + y/3)M, eq. (|3.93| ) and eq.( 3.94 ) have a solution simultaneously. In this case, 
however, the solution does not correspond to the absolute minimum of the potential V and the potential at 
the true minimum becomes negative. Therefore, we do not consider this case. 
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By solving eg. ( ^.931 ), vacuum-expectation value of the Polonyi field 0p is given by 



)= V3-1 M, 



(3.95) 



and one can obtain the vacuum-expectation value of the superpotential Wp and the order 
parameter Fp; 



{Wp) = fi^M, 
{Fp) = VSe^-^/i^. 



(3.96) 
(3.97) 



From the vacuum-expectation value of Wp, the gravitino acquires a non- vanishing mass 

,2 



fns/2 



,2-^3 /i 



(3.98) 



Notice that for the given gravitino mass ^3/2, the scale of the condensation of Fp (which is 
the same order of /i^) is determined; {Fp) ~ 0(m3/2M). 

The mass eigenvalues m^p^ and m^p^ for the scalar fields are obtained by expanding 
Polonyi field 0p around its vacuum-expectation value; 



(,p = (y3-l)M+-^(0Pi + 



t(pP2) . 



(3.99) 



Substitute this to the potential ( |3.92|) , the masses are given by 



^Ipi = 2V3"^3/2. ^1 



(4 - 2V3) 



m 



3/2 5 



(3.100) 



while the superpartner of 0p is a goldstino and absorbed in the gravitino. 

Now, let us consider the coupling between hidden and observable sectors. The hidden sec- 
tor couples to the observable sector if we introduce a superpotential Wohs for the observable 
sector; 



w = Wp + w,u<Pobs) 



(3.101) 



where we denote the particles in the observable sector by 0obs- Assuming the minimal Kahler 
potential for (pohs, the scalar potential is given by 



exp 



M2 



^W + 



dWp 



+ E 



'ohsyy^dWoh^ 



M2 



obs 



In the potential ( |3.102D , all interactions between the hidden and the observable sectors are 
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suppressed by powers of M^^. 

In order to derive a low energy effective tlieory for tlie observable sector, we take a limit 
M — i> oo with gravitino mass ^3/2 = {e^^'^^^W/M'^) fixed. (This is sometimes called flat 
limit.) Then, one can obtain a potential for the observable sector; 



bbs I 



E 



obs 



obs 



+ ^1/2 E 



^obs 



+m3/2 



where 



E '/'obsS^ + {b*{a + 6) - 3} lyobs + h.c. 



^obs 



(3.103) 



W^obs = (e^/^^^^)iyobs, 
'{dWp/d(j)p)M\ 
/ 




(3.104) 

a- = { ^T^- — ) = 1, (3.105) 

(3.106) 

The above potential is nothing but a potential for the softly broken (global) SUSY. The first 
term in eg. ( |3.10^) is a scalar potential in a global supersymmetric model, while the rest can 
be regarded as soft SUSY breaking terms for scalar bosons. Especially, in this model all the 
scalars in the observable sector receive an universal SUSY breaking mass which is fixed by 
gravitino mass ^3/2- 

In order to keep the masses of squarks and sleptons at O(100GeV — ITeV) so that the 
hierarchy problem is solved, the gravitino mass m3/2 is taken smaller than O(lTeV). From 
this phenomenological requirement, (Fp) is determined; 

(Fp) ~ (10^^ - 10^^)GeV, (3.107) 

provided 77^3/2 ~ O(100GeV - ITeV). 

Finally we will comment on gauge fermion mass. In a model with the minimal ki- 
netic function for gauge multiplet; f(ab) = ^(ab), gauge fermion remains massless. To gener- 
ate a non-vanishing gauge fermion mass term, we assume a non-minimal kinetic function; 
idf(ab}/d4>p) 7^ 0. Then with a SUSY breaking in the hidden sector, gauge fermion acquires 
mass from the interaction in the total supergravity lagrangian ( p.32[ ); 



C 



XX 



df(ab) 



X^^h^''^ + h.c. 



(3.108) 



The simplest example for the non-minimal kinetic function is 



(3.109) 
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where k is a dimensionless coupling parameter. Then, the gauge fermion mass is given by 
K (Fp) /2M, which is of the order of the gravitino mass for ^ ~ ^(1)- 



3.7 Mass of a scalar field in the hidden sector 

As we have seen in the previous section, masses of the Polonyi fields 0pi and (j)p2 are 
of the order of the gravitino mass ^3/2- In fact, this is not an accidental case, i.e. in 
supergravity there exists a scalar field with mass of order m3/2 in a wide class of models with 
the following features. 

• At the stationary point of the scalar potential V, cosmological the constant vanishes. 

• SUSY is broken by a condensation of a F-term in the hidden sector. (In this case, 
the vacuum-expectation value of (F) is 0{m^i2M) so that the cosmological constant 
vanishes.) 

• Cut-off scale of the model is of the order of the gravitational scale M. 

As we will see below, a scalar boson with mass of order 7713/2 exists in a model with these 
conditions. 

The existence of such a scalar field can be seen by using only simple order estimations. 
Assuming the cosmological constant to vanish, the condensation of F-terms are constrained 



{F') < 0{{e^'^^^^^W/M)) ~ 0(m3/2M) (for all i). (3.110) 
Combining this condition with the stability condition (dV/dcp) = 0, we get p3, |2P 



Me^/^^'GijF^) < 0{ml/2M), (3.111) 

where we have used that the cut-off scale of the model is of order M, i.e. {Gijk*) ^ M~^. 

The mass squared matrix of the scalar fields are given by the second derivative of the 
scalar potential V. Especially, the diagonal part is given by 



d(j)^d(j>> 



- {{Me''/'^''G,k){g''*K,'rm){Fn + h.c.) 

+0{ml/2). (3.112) 



^In order to obtain a correct normalization of kinetic terms, we take {Gij*) — Sij* . 
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Denoting the chiral multiplet which breaks SUSY {(p^ , , F^) , the vacuum expectation 
value of is given by (F-^) ~ 0(m3/2M). Then, constraint (|3.111| ) leads to 



(Me«/2*^'Gx.)<0(m3/2), (3.113) 

and hence 

<Oiml,). (3.114) 



d(f)^d(j) 



X* 



From the above, we can conclude that the mass-squared matrix of scalar fields has a eigen- 
value of order m\i2 or less. We should note here that such a scalar field may cause serious 
cosmological difficulty |^0| (so-called Polonyi problem) which will be discussed in Chapter |[ 
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Chapter 4 

Feynman rules for the gravitino 



In the previous chapter, we have seen that the gravitino, which is the gauge field associ- 
ated with local SUSY invariance, plays a crucial role in supergravity. In order to see physics 
of the gravitino, we must discuss a field theory for it. 

In supergravity, however, there exist non-renormalizable interactions which are sup- 
pressed by (or its higher power), and hence we cannot calculate loop effects by us- 
ing the full supergravity lagrangian given in the previous chapter. Furthermore, because of 
these non-renormalizable interactions, the Born-unitarity breaks down at high energy scales 
of order M. Therefore, we regard the supergravity as a low energy effective theory which is 
appropriate for energy scales much below the gravitational one, and only consider tree level 
processes. 

In this chapter, we derive Feynman rules for the massive gravitino field by using the 
supergravity lagrangian ( |3.32| ). For this purpose, we first solve the field equation for massive 
Rarita-Schwinger field, and then we discuss the quantization of the free gravitino field. 
Hereafter, we only consider the nature of the gravitino in (nearly) fiat space-time, and hence 
we restrict the background metric in that of fiat space-time g^^ = diag(l, —1, —1, —1). 



4.1 Four-component notation for fermions 

For our following arguments, it is more convenient to use a four-component spinor for a 
fermionic field rather than two-component one used in the previous chapter. Therefore, we 
briefiy introduce it before quantizing the gravitino field. (For our notations and conventions, 
see Appendix 

The four-component spinors ip (in chiral representation) can be constructed from two 
component spinors C, and rj; 

^-(h)^ (4.1) 
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Corresponding to this, 7-matrix consists of 2 x 2 matrices a'^ and a'^; 

In the above representation, we define four-component spinors for the gravitino, 
gaugino, A^^-*, and chiral fermion, Xh'^ 

. *r-(*V.*.). (4.3) 

>''■'"'' ^(j) ■ a'"'^(A.a), (4.4) 

xf'^(l) . x'r'^m). (4.5) 

Notice that iljj!^^ and A*.*^-* satisfies the Majorana condition 



where C is the charge conjugation matrix (see Appendix ^), while Xr is chiral fermion 
with a definite chirality; 

^(l + 75)xS?^ = xi?\ ^(l-75)xi?^ = 0, (4.7) 
where 75 = i7''7^7^7'^. 

4.2 Wave function for massive gravitino 

We will start by discussing the wave function for massive gravitino. The lagrangian 
for the free gravitino field is given in the full supergravity lagrangian (|3.32|) . In the four- 
component notation, it can be written as 



Cmrs = -^e'^^^^^^757.5p^. - ^"^3/2^^ [7^ Y] (4- 



where we have dropped the index (M) for the gravitino field, for simplicity. (Hereafter, we 
use only the four- component notation otherwise mentioned, and drop indices (M) and (D) 
in order to avoid complications due to too many indices.) By using the Majorana condition 

—T „ 

for the gravitino field; if)^ = Cip , the lagrangian (|4.8| ) becomes 



Cmrs = \e^'""'V^C^i^i.d,^, + ^m3/2<Ct [7^ 7^ (4.9) 



This lagrangian is our starting point 
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Varying the above lagrangian for t/'^, we can obtain the field equation for the free gravitino 



field; 



d 



MRS 



From this, the following two equations are derived; 



(4.10) 



(4.11) 
(4.12) 



Notice that the former (latter) equation can be obtained by operating 5^ (7a7^i) on eq. (|4.10|) . 
For a later convenience, we derive one more equation by multiplying eq. ( [4.12|) by 7^; 



t {jdYi^^ - 7VV'a) + 3m3/27'^V'M = 0. 



(4.13) 



For the massive gravitino (m3/2 7^ 0), eq.( |4.11| ) - eq.( |4.13| ) yields the following simple equa- 
tions; 



0, 
0, 
0. 



(4.14) 
(4.15) 
(4.16) 



The solutions to the above equations can be constructed by using the wave function u for 
spin i field and the polarization vector for spin 1 field. In solving the field equations for 
the gravitino field ip^, it is convenient to consider the wave function in the momentum 
space, '^p^i ~ e'^^^ip^. Then, the solution to eq. ( [4.14| ) - eq. ([4.16|) is given by ||31| 



^.(P,A)=E((^'£)(1'^) 



A) )n(p, s)e^(p,m). 



(4.17) 



where {{■^, ^){l,m) \ (|,A)) is the Clebsch-Gordan coefficient, whose value is shown in Ta- 
ble O. 



The wave function for spin | field u{p, s) is the solution to the ordinary Dirac equation 
(in momentum space) with a definite helicity (s = ±1); 



7713/2) m(p,s) = 0, 
(nS)M(p, s) = su{p,s) 



(4.18) 
(4.19) 



where nS is the helicity operator (see Appendix The explicit form of u{p, s) in the Dirac 
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m = —1 


m = 


m = +1 


s = -l 


1 


^2/3 


\/l/3 


s = +1 


\A73 


^2/3 


1 



Table 4.1: Clebsch-Gordan coefficients for the case A = | + m; ((|, |)(l,m) 
Notice that the coefficient with A 7^ | + m vanishes. 



;|,f + m)). 



representation is given in Appendix 0, and we take the following normalization condition 
on m(p, s); 



u{p,s)u{p,s ) = 2m3/2(5^,'. 



(4.20) 



For the momentum vector = {E, |p| sin ^ cos 0, |p| sin ^ sin 0, |p| cos^) of a massive 



particle {p^p^ = m\/2 7^ 0), the polarization vectors take the following forms; 



e/.(p,l) 
eM(P>0) 



1 



"^3/2 
-1 



(0, cos6' COS0 — i sin0, cos6'sin0 + icos0, — sini 



|p|, — E'sin^cos^, — E'sin^sin^, —Ecosi 



e^(p, — 1) = ^p(0, cos6' COS0 + i sin0, cos 6' sin — i cos ( 
v2 



sm I 



(4.21) 
(4.22) 
(4.23) 



which are normalized as 

e;(p,m)e'^(p,m') = (4.24) 
Notice that the polarization vectors ( [4. 211 ) - ( |4.23| ) satisfy the following condition; 



p''e^(p, m) = p'^e* (p, m) = 0. 



(4.25) 



By using eq. (^4.18|) , eq. (|4.19|) and eq. (|4.25| ), one can easily check that defined in 
eq.( |4.17| ) obeys the following equations; 



Y^,{pA) = 0, 

p''^^(p,A) = 0, 
(/»-m3/2)^^(p,A) = 0, 



(4.26) 
(4.27) 
(4.28) 



and hence the wave function in the coordinate space ip^ ^ e ^^^ip^ satisfies eq.( [4.14 ) - 
eq.(|4l6|). 

The explicit form of '0'^(p, A) depends on representations of the 7-matrix, and we do 
not write it down because it is complicated but almost useless. Instead of that, we give 
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some useful identities for ipfj,{p, A), which hold irrespective of the representation of 7-matrix. 
Normalization of i'^ip, A) is fixed by eg. (|4.20|) and eg. (|4.24| ); 



(4.29) 



Furthermore, one can obtain the following identity which will be useful in deriving the 
momentum operator given in eg. ([4 .541) for the gravitino field; 



V^^(p,A)7i.V^^(p,A') = -2pJ^^'. 
For the helicity sum of the gravitino field, the following formula exists; 

PfiAp) = Z!^''(P'^)^/.(P'^) 

A 



m3/2 



9^lv 



mt 



3 \ ^^" 



PuPX 1 a A 
9uX — 7 7 

mi 



''3/2 / '-' \ '"'3/2 / \ "^3/2 

P^uip) obeys the following eguations which correspond to eg.( [4.26|) - eg.( [4.28|) 



(4.30) 



(4.31) 



7'^P^,(p)=P^,(p)7^ = 0, 
p'-P.M = PMp" = 0, 

- "^3/2) Pf^uip) = Pf^uip) ii) - "^3/2) = 0. 



(4.32) 
(4.33) 
(4.34) 



General solution to eg.( [4.14| ) - eg.( |4.16| ) can be expanded by the mode function given in 
eg.( |4.17| ) and its charge conjugation; 



^0^^ E {e^^>.(P, A)apA(t) + e-P>J(p, A)at.,(t)} 



(4.35) 



where po = y |Pp + "^3/2' ^^'^ '^pa(^) is the expansion coefficient whose time dependence is 
determined by the field eguation ( ^4.16| ); a-px{t) ~ e"*^"*. Notice that because of the Majorana 
nature of the gravitino field tp^^ the coefficient of the ?/'^(p, A)-term in eg.( [4.35| ) is not an 
independent variable, but it is to be a}px{t). One can easily check that the gravitino field ijj^ 
in eg.( f4.35| ) satisfies the Majorana condition; ifj^j^ = ip'^. 



4.3 Quantization of free massive gravitino field 

In order to obtain Feynman rules for the gravitino, we have to guantize first the free 
gravitino field. In this section, we discuss the canonical guantization of the free gravitino 
field. Since in this thesis, we only consider spontaneously broken local SUSY models in 
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which the gravitino acquires a non-vanishing mass term, we only deal with the case of the 
massive gravitino. 

Constraints on the physical mode of the gravitino field are obtained in the previous 
section, and they are given in eg. ( [4:. 14 ) and ( [1.15| ). The mode expansion of the physical 



degrees of freedom is given in eq.( [4.35| ), and as in the usual quantization procedure, we 
regard the coefficients apA(t)'s and apA('^)'s as dynamical variables and derive commutation 
relations among them. 

As a first step in the quantization procedure, we derive canonical momenta. Time deriva- 
tives of apx{t) and cipxit), which we denote apx{t) and ap;^(t), are present in the lagrangian 
as 



MRS = / " -^-^MRS 
3, 



1 . f d^p 



2 J (27r)32po X 
+ (terms without dp^, ctpA)- (4.36) 



Differentiating lagrangian ( |4.36 ) with respect to dpx(t) or d^-^it), one can obtain canonical 
momenta IIpA and IIpA; 

HpA = -Y- = ^o /o ^3o %aW, 4.37 

bapx 2 (27r)'^2po ^ 

Hda = 1 — = — cidaU^)- 4.38) 

bd\,^ 2 (27r)^2po "^^^ ^ ^ 

Since dp\{i) and a'px^) cannot be expressed as functions of canonical momenta, eq.( |4.371 ) 
and eq. ([4.38|) are regarded as primary constraints on this system; 

For this system, Poisson bracket {■ ■ - jp is defined as in the usual case; 

bF\ / 6G \ f 6f\ / 6G \ 



{F,G}p ^ Jd'pY: 



R VOiipA/ L yapX j R \«-^^pA/ L 

where the index L (R) represents left (right) derivative. Especially for the dynamical vari- 
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ables apx{t), cipxi^) canonical momenta IIpA, HpA, Poisson brackets are given as 

{apA(t),np,v}p = 5,v5(p-p'), (4.42) 
{at,,(t),np,v}p = 5,y5(p-p'), (4.43) 

and those for constraints $pA and $pA can be obtained as 

^ f '5aa"^(p-p: 



(27r)32po I 5aa'^(p-p) 



(4.44) 



Notice that the above matrix C (^p, A;p',A'^ has its inverse, and hence time evolutions of 
the primary constraints ( ^4.39| ) and ( [4.40| ) do not induce secondary ones. 



In order to quantize this constrained system, we introduce Dirac bracket {■ ■ ■}■£,; 
{F,G}^ ^ {F,G}^- J d'pd'p Y.{F,Kx}pC-' {p,X-p\x') g}^ 

XX' 

= {F, G}p + {271 f t [d'pYl 2Po {F, %x}p {$pA, G} 

X 

+ {2n ft I f/3p^2po{F,$pA}p{$pA,G}p, (4.45) 

where $pA represents both $pA and $pA, and F and G arbitrary variables. Then, this system 
is quantized by replacing the Dirac bracket by (anti-)commutator; 

t {F, G}d [F, G} = FG- (-l)l^ll^lG'F, (4.46) 

where \F\ = 1 if F is Grassmann-odd and \F\ = if F is Grassmann-even. Especially, 
commutation- relations among apA(t) and a.pA(^) are obtained as 

{apA(t),a|,,,,(t)} = (27r)'2po5AA'<^(p-p'), (4-47) 
{apA(t),ap,y(t)} = {al,{t),al,^,{t)} = 0. (4.48) 

As in the ordinary procedure, we can get a hamiltonian for this system; 

H = Jd^pY. (HpAOpA + TTpAaJ^A) - / c^^xA 



-MRS 



/ 1— 1 — 

d^p 

W)Wo A 



/74l^E^oaiA(t)apA(t). (4.49) 
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By using this hamiltonian, equations of motion for ap\(t) and a^pxi^) derived; 



and these equations can be easily solved; 

where we denote apx{t = 0) and a^xit = 0) as apx and a^x- 

Then, by using apx and a^x, the field operator can be expanded as 



(4.50) 
(4.51) 

(4.52) 



(27r)32po , 

Furthermore, the momentum operator for the gravitino field is given by 



(4.53) 



Pn 



(27r)^2po 



(4.54) 



Fock space for the gravitino field is constructed by operating the creation operator a^x 
on vacuum state |0) which satisfies the condition apA|0) = 0. Especially, one particle state 
of the gravitino with momentum p and helicity A is defined as 



|P,A)=a;,|0). 

As one can see, this state satisfies the following relations; 

{p,X\pA') = i2nf2po6xx'S{p-p), 
Pf^lpA) = Pt^\p,^), 



(4.55) 



(4.56) 
(4.57) 



where is the momentum operator defined in eq.( |4.54| ). Notice that the normalization 
condition ( [4.56|) on the one particle state ( [4.55| ) is Lorentz invariant since PoS{p — p) is a 
Lorentz invariant variable. 

From the above arguments, we obtain Feynman rules for the external gravitinos; for the 
gravitino with momentum p and helicity A in initial the state, we should assign ip^ip, A) or 
ip*^' ^{p, A), and that in the final state, 'ip^ip, A) or ip'^{p, A). 
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4.4 Interactions of the gravitino 



In the previous section, we have quantized the massive free gravitino field, and obtained 
Feynman rules for the gravitino in external line. Our next purpose is to discuss the Feynman 
rules for the interactions of the gravitino field. Many interaction terms concerning the 
gravitino field ip^ exist in the full supergravity lagrangian (|3.32|) , but most of them are 



irrelevant for our present purpose. This fact arises from the following two reasons. 

• Since we are considering the processes with the energy scale much below the gravita- 
tional one; a/s <^ M, contributions of the Feynman diagrams with higher dimensional 
operators are suppressed by factor ~ \fs/M ^ 1 (or its higher power). As we will see 
later, the dimension of relevant operators is always five, and interaction terms with the 
dimension higher than six are not important for us. 

• In our analysis in Chapter |^ and Chapter 0, gravitinos only appear at external lines in 
the Feynman diagram. Then, interaction terms which contain 7^'?/'^ or ip^"^^ are not 
necessary since they vanish due to eq. ([4.26|) . 

In the following arguments, we take only the relevant interaction terms into account and 
ignore contributions from other terms. 

The most important interaction terms come from the coupling between the gravitino field 
and the supercurrent, which is given in eq. (|3.65|) . In the four-component notation (and in 
the fiat space-time), these terms are written as 

fclT-.TlT-AWFiJ. (4.58) 



I 



8M 

From this, we construct Feynman rules for the interaction of the gravitino field ipf^ with 
matter fields 0, ^^^id A. In Fig. [4.1| , we show the Feynman rules derived from the 



lagrangian ( 4.58 ). In the following analysis, we use these rules and ignore other higher 
dimensional interactions. 

As a simple example of applications of these Feynman rules, we calculate the decay 
rates of the gravitino. The dominant decay modes of the gravitino are ipfj_ —>■ X + Afj_ and 
0* + X* (and — i> 0** + x') if these processes are kinematically allowed. Feynman 
diagrams for these processes are shown in Fig. [4.2| . 

For the process ^/^^j ^ A + A^, the invariant amplitude can be obtained as 

M{^, ^X + A,) = ^q^eA {g^'Y - g'^Y) (4.59) 
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By using the mode sum for ip^ given in eq.( |4.31| ), eg. ( [4. 591 ) becomes 



1 



Asm 
4 

_ "^3/2 
2M2 




(4.60) 



where m\ is the gauge fermion mass, and we have assumed that the gauge field is massless. 
Notice that in deriving eq.( |4.60D , we have taken an average of the hehcity of initial gravitino. 
By using eq. (|4.6CI| ), the decay rate is given by 




= 3^ll{^-|£^l Mi + ^(£^n. (4.61) 

where P/ is the three momentum of the particle in final state. Notice that if the gauge group 
G is non-abelian, the above decay rate should be multiplied by the rank of G in order to 
calculate a total decay rate. Especially in the case ^ mx, the gravitino lifetime is 
approximately given by 

Tsd^, ^ A + A^) ^ 4.0 X 10' X (y^^) ' sec. (4.62) 
For the process ipij. ^ 4>^ + X^, the invariant amplitude is given by 

MiiJ, ^ f + xi = ^^X7^/(l + 75)^M, (4.63) 
and by using a similar method to in the previous case, the decay rate is obtained as 

where is the mass of (p field, and we have assumed that x is a massless fermion. Notice 
that the decay rate for the process ipfi 0'* + X* is equal to the one given in eq. ( 



Another application is to calculate the gravitino production cross sections. As one can 
easily see, the most important processes are those of one gravitino production since the 
vertices having a gravitino are suppressed by (or its higher power). Amplitudes for 

the gravitino production processes are obtained by combining the Feynman rules given in 
Fig. [4.1| with those derived from the ordinary (global) SUSY lagrangian. We have calculated 
the total cross section for the dominant processes of helicity ±| gravitino production, and 
the results are shown in Table W^- 
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Process 



a 



(^V647rM2)x 



(A) 
(B) 
(C) 

(D) 
(E) 
(F) 
(G) 
(H) 
(I) 
(J) 



A" + - 

+ A(^) 
A"^ + 0, - 

Xi + 0* - 
A{«) + 

Xi + Xj - 



+ A"" 
^i^ + Xj 

> ^ + A" 

^ + A(^) 
^i^ + Xj 



abc 



/3) / 

fabc 2 {_(3/2) + 2 ln(2/5) + 5 - (1/8)52} 



ji 
cabc 



pbc {-(62/3) + 16 ln[(2 - 5)/ 5] + 225 - 25^ + (2/3)5^} 



''pa 
ji 
2 

ji 



(8/3) 
(16/3) 



{-2 + 21n(2/5) + 5} 
{-6 + 81n(2/5) +45- 

2 

ji 

2 



;i/2)52} 



ji 



Table 4.2: Total cross sections for the helicity ±| gravitino production processes. Spins of 
the initial states are averaged and those of the final states are summed. /"'"^ and T^j represent 
the structure constants and the generators of the gauge groups, respectively. Notice that for 
the processes (B), (F), (G) and (H), we cut off the infrared singularities due to the t-, u- 
channel exchanges of gauge bosons, taking a small but a positive parameter 5 = (1 ±cos 9)min 
where 6 is the scattering angle in the center-of-mass frame. 



4.5 Effective lagrangian for light gravitino 

In spontaneously broken SUSY models, the massless gravitino field '^p^l acquires mass by 
absorbing goldstino modes. Before becoming massive, the gravitino only possesses helicity 
±1 modes, and the goldstino provides helicity ±| modes of the massive gravitino field. This 
fact suggests that the helicity ±| mode of the gravitino field behaves like a goldstino. In 
fact, if the gravitino mass m3/2 is much smaller than the mass differences between bosons 
and fermions in the chiral and gauge multiplets, the above argument is valid and we can 
obtain an effective lagrangian for the relativistic gravitino field of helicity ±| components. 
In this section, we will derive the effective lagrangian for the light gravitino field and 
apply it for calculating some processes. 

For the case y/s ^ ^3/2, the wave function of the gravitino of helicity ±| components 
is approximately proportional to p^/m^i2 where is a momentum of the gravitino. In this 
case, the hehcity ±| component of the gravitino field can be written as 

V 3 7713/2 

where tp represents the spin ^ fermionic field which can be interpreted as the goldstino. 
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Substituting eq.( [4.65| ) into the gravitino interaction lagrangian ( [4. 581 ), we obtain the effective 
interaction lagrangian for the goldstino components ip- 



Using the replacement ([4.65 ), the gravitino interaction lagrangian ( [4.58 ) becomes 



eff 



4x/6Mm3/2 ^ ^ ' ' ' J ' " 



+ (total derivative), (4.66) 

where we have used the relation jdtpfj, = since we are considering processes with a light 
gravitino in the external line. 

The terms in the right-hand side of eq.( [4.66[ ) are dimension six operators and one would 



expect that the processes involving helicity ±i gravitinos have bad high energy behaviors. 
In supergravity, however, this is not the case since the leading high energy behavior cancels 
out in the total amplitude. For example, for the helicity ±^ gravitino production process 

are s-, t- and w-channel diagrams (see Fig. [4.3| ), whose leading 

terms are given by 



^- - H/|;^(r"),<V'/A- (4.67) 



where p' is the momentum of 0*. As one can see, they cancel out in the total amplitude. This 
can be understood in the following way. The helicity ±| components of the gravitino field is 
the unphysical in the symmetric phase, and hence total amplitudes with helicity ±| gravitino 
at external line should vanish unless SUSY is broken. That is, the total amplitude for the 
helicity ±| gravitino production should be proportional to some SUSY breaking parameters. 
Thus, the leading high energy behavior, which is independent of SUSY breaking parameters, 
cancels out in the total amplitude. 

From this fact, we can obtain an effective lagrangian for the helicity ±| light gravitino 
field by replacing all the derivatives operated on 0, x, and A in eq.( [4.66| ) by the masses m^, 
m^, and nix of the corresponding fields. When these derivatives are operated on external 
lines, this replacement is easily justified. On the other hand, if they are operated on prop- 
agators, they become internal momenta and seem to make the high energy behavior worse. 
As mentioned above, however, this bad high energy behavior cancels out, and hence the 
leading high energy behaviors can be subtracted from the propagators with derivatives in 
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Process 




{-12 + 161n(2/5) + 85-^2} 



{-22 + 241n{(2 - 5)/ 5} + 245 - 35^ + 
{-2 + 21n(2/5) + 5} 
^ {-28 + 32 ln(2/(5) + 165 - 5'^} 



Table 4.3: Total cross sections for the helicity ±| gravitino production. Spins of the initial 
states are averaged and those of the final states are summed, f"'^'^ and T""- represent the 
structure constants and the generators of the gauge groups, respectively. Notice that for 
the processes (B), (F), (G) and (H), we cut off the infrared singularities due to the t-, u- 
channel exchange of gauge bosons, taking small but positive 5 = (1 ± cos 6')mm where 6 is 
the scattering angle in the center-of-mass frame. 
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(4.70) 
(4.71) 
(4.72) 



Thus, the derivatives in the interaction lagrangian ( [4. 661) can be replaced by the appropriate 
masses related to the SUSY breaking, and the lagrangian ( [4. 66]) becomes 



eflf 



i{m1 -m^) /_ 



1/37713/2 M 



-tmx — 



8VQrn3/2M 



(4.73) 



As one can see, in the SUSY limit [i.e. — — > and mx — > 0), the above lagrangian 
vanishes and the helicity ±| modes of the gravitino field decouple from the theory. 

In the case of a light gravitino, it is more convenient to use the effective lagrangian ( [4.73]) 
rather than the full lagrangian ( |4.58| ) . Feynman rules derived from the lagrangian ( [4.73| ) are 
shown in Fig. W^. 
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By using the effective lagrangian ( [4. 731 ) , we have calculated decay rates and cross sections 
for some processes including a light gravitino. The decay rate of the process X ^ ip + is 
given by 



2 



r(A-^. + ^,) = -^^ l-(^yK (4.74) 



and that of the process (p ^ ip + 




^(^-* + ^-'-4S^^l-|-^n. (4.75) 

Furthermore, total cross sections of some light gravitino production processes are calculated 
and the results are shown in Table Notice that the interaction of the helicity ±^ gravitino 



becomes stronger as the gravitino mass ^3/2 becomes lighter, as one can see in eq.( |4.74| ), 
eq.( |4.75| ) and the total cross sections given in Table 



The results obtained above can be used when the gravitino mass 7713/2 is much smaller 
than the SUSY mass splitting in observable sector. Such a situation seems to conflict with 
the super-trace formula obtained in the previous section. But the super-trace formula is a 
tree level relation, and the mass splitting in the observable sector may become larger than 
the gravitino mass due to various radiative corrections. In fact, in some models based on 
no-scale type Kahler potential |]10|, the gravitino mass may become much smaller than the 



electroweak scale. In any case, if the gravitino mass is smaller than the mass splitting in 
other multiplets, we can use the results obtained in this section. 
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Figure 4.1: Feynman rules for the interactions of gravitino. 
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4 

II 




(a) 



(b) 



Figure 4.2: Feynman diagrams for the processes (a) ■^^ — > A + A^, and (b) ■0// 0' + 





Figure 4.4: Feynman rules for the interactions of a light gravitino. 
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Chapter 5 

Phenomenology of the gravitino : 
overview 

Having Feynman rules with the gravitino field in the previous chapter, we are now at 
the point to discuss phenomenology of the gravitino. As we will see below, the gravitino 
is almost nothing to do with collider experiments because its interactions are extremely 
weak. But if we think of cosmology, effects of the gravitino may become very significant. 
Our main purpose is to investigate the effects of the gravitino on the inflationary universe 
quantitatively. Before doing this, we will survey the phenomenological implications of the 
massive gravitino. 

5.1 Collider experiments with the gravitino 

In this section, we will give brief comments on collider experiments with the massive 
gravitino. If the gravitino mass is comparable to the masses of SUSY particles {i.e. squarks, 
sleptons and gauginos), gravitino production cross sections are extremely small since the 
interaction terms for the gravitino are suppressed by powers of M~^. For example, the 
helicity ±| gravitino production cross section in e+e" collider experiments is estimated to 
be 

a(e+ + e- ^ V^, + A) = + ^gj) ~ 1 x IQ-^W, (5.1) 

where we have assumed that ^/s is much larger than the gravitino and gaugino masses. 
Compared with the luminosity of LEP ~ lO^^cm^^sec"^, or even with a design luminosity 
of JLC ~ lO^^cm^^sec"^, we have no hope to have a signal of gravitino production events. 

In the case of light gravitino, however, helicity ±| modes of the gravitino interact strongly, 
and hence one may expect that some signals of the gravitino may be detectable if the gravitino 
mass is sufficiently small. Bound on the light gravitino mass from collider experiments 
depends on the mass spectrum of superparticles. The most stringent bound on the gravitino 
mass is derived for the case where there exists a neutralino lighter than Z°-boson. In this 
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case, Z°-decay may contain a single photon balanced by missing transverse momentum in 
the opposite hemisphere [Q. Such a single photon is expected to arise from the process 



^ i'fi + X: X ^ i'fi + where x is the neutralino. Both of the decay rates for these 
processes are proportional to rn^^^, and hence such a single photon event is sizable if the 
gravitino mass is extremely small. This process is analyzed in ref. |Q, and the gravitino 



mass smaller than ~ 10^^'^GeV is excluded by the LEP experiments. Notice that this bound 
(^3/2 ~ 10^^'^GeV) is comparable to that obtained from cosmological considerations. If the 
neutralino mass is larger than the Z'^-boson mass mz, we cannot use the above arguments, 
and less stringent bound may be derived. 

5.2 Cosmology 

Contrary to collider experiments, the mass of gravitino is severely constrained if we 
assume the standard big-bang cosmology.Q If the gravitino is unstable, it may decay after 
the big-bang nucleosynthesis (BBN) and produces an unacceptable amount of entropy, which 
conflicts with the predictions of BBN. In order to keep the success of BBN, the gravitino 
mass should be larger than ~ lOTeV as Weinberg first pointed out [ill]. In a model with the 
gravitino mass larger than ~ lOTeV, the gravitino decay processes produce a large amount 
of entropy and dilutes the baryon number density of the universe. This requires that the 
baryon-to-photon ratio before the decay of the gravitino should be extremely large so that 
the present value of the baryon-to-photon ratio is given by O(10~^ - 10~^°). Furthermore, in 
a model with i?-parity invariance, the gravitino decay also produces an unacceptable amount 



of LSP which conflicts with the observations of the present mass density of the universe |]12 
(These problems are sometimes called "gravitino problem".) 

Meanwhile, in the case of stable gravitino, its mass should be smaller than ~ IkeV not 



to overdose the universe [|T3[. Therefore, the gravitino mass between ~ IkeV and ~ lOTeV 



conflicts with the standard big-bang cosmology. 



However, if the universe went through inflation we may avoid the above con- 



straints since the initial abundance of the gravitino is diluted by the exponential ex- 



pansion of the universe. But even if the initial gravitinos are diluted, the above problems 
still potentially exist since gravitinos are reproduced by scattering processes off the thermal 
radiation after the universe has been reheated. The number density of the secondary grav- 
itino is proportional to the reheating temperature and hence upperbound on the reheating 
temperature should be imposed not to overproduce gravitinos. Therefore, even assuming 
the inflation, a detailed analysis must be made to obtain the upperbound on the reheating 
temperature. 



^In the "standard cosmology", we adopt the following two points. The first point is that the light 
nuclei were synthesized through the (almost) standard big-bang nucleosynthesis, and the second is that the 
density parameter SI of the present universe is smaller than 0(1). (For the big-bang nucleosynthesis, see 
Appendix ^.) 
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For the case of the unstable gravitino, the cosmological considerations on regeneration 
and decay of the gravitino has been done in many articles 0, ^ |3^, ^ . These 



previous works show that the most stringent upperbound on the reheating temperature comes 
from the photo-dissociation of light nuclei (D, T, ^He, "^He). Once gravitinos are produced in 
the early universe, most of them decay after BBN since the lifetime of the gravitino of mass 
O(100GeV — lOTeV) is O((10^ — 10^)sec). If the gravitinos decay radiatively, emitted high 
energy photons induce cascade processes and affect the result of BBN. Not to change the 
abundance of light nuclei, we must have a restriction on the number density of gravitinos, 
and this constraint is translated into the upperbound on the reheating temperature. 

In order to analyze the photo-dissociation processes, we must calculate the following two 
quantities precisely; the number density of the gravitinos produced after the reheating, and 
the high energy photon spectrum induced by radiative decay of the gravitinos. But the 
previous estimations of these values are incomplete. As for the number density of gravitino, 
most of the previous works follow the result of ref . |^ , where the number density is under- 
estimated by factor ~4. Furthermore, in many articles the spectrum of high energy photon, 
which determines the photo-dissociation rates of light elements, is calculated by using a sim- 
ple fitting formula. In this thesis, we calculate these two quantities precisely and find more 
precise upperbound on the reheating temperature. 

If the gravitino is stable, the situation changes. If a very light gravitino of mass m^i2 ~ 
0(1 keV) was once thermalized, the critical density of the universe is easily obtained, and 
hence it constitutes a dark matter. For a heavier gravitino, one has to consider inflation. 
We will find that the requirement that the gravitinos produced after the inflation should 
not overdose the universe (the closure limit) places a severe upperbound on the reheating 
temperature Another stringent constraint comes from the photo-dissociations of light 
nuclei by high energy photons produced through radiative decays of the next-to-the-lightest 
superparticle into the gravitinos, which excludes the certain range of the gravitino mass. 

In the following chapters, we will investigate these effects quantitatively. 



5.3 The Polony i problem 

As we have seen in Chapter |], there exists a scalar boson with mass of the order of the 
gravitino mass m3/2 in a wide class of models based on supergravity. (Below, we call this 
scalar field "Polonyi field" 0p and denote its mass m^p.) If one assumes that the gravitino 
mass is of the order of the electroweak scale from the naturalness point of view, it has been 
pointed out that the Polonyi field 0p causes serious cosmological difficulties (so-called 
"Polonyi problem") as gravitino does. Though this is not directly related to our present 
purpose, we briefly discuss this problem because of its significance. 

The Polonyi problem stems from the fact that the Polonyi field takes its amplitude of 
order M at the end of inflation, i.e. there is a Bose condensation. Such a condensation cannot 
be eliminated by inflation since the equation of motion for the Polonyi field is coupled to 
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that of the inflaton field. Thus, in general the potential minimum of the Polonyi field at zero 
temperature is not a stationary point during infiation, and hence the Polonyi field does not 
sit at the origin. Furthermore, it should be noted that scalar fields deviate from the origin 
due to the quantum fiuctuation during the infiation. 

As the temperature drops and the expansion rate H of the universe becomes compara- 
ble to m^p, the Polonyi field starts its oscillation, and subsequently dominates the energy 
density of the universe until it decays. The decay rate of the Polonyi field is estimated 
as 0{m^p/Mpi) ~ 0{my2/^pi): ^ind hence it may cause cosmological difficulties as in the 
gravitino case. That is, the Polonyi field decay releases tremendous amount of entropy and 
dilutes primordial baryon asymmetry much below what is observed today (which is some- 
times called the entropy crisis). Furthermore, the decay of the Polonyi field destroys the 
successful scenario of the BBN if it occurs after the BBN starts. In fact, the Polonyi prob- 
lem is more serious than the gravitino problem, since it cannot be solved by infiation. As 
a result, a wide class of models with the SUSY breaking in hidden sector are excluded by 
cosmological arguments. 

It has been pointed out |^ that the first problem [i.e. the entropy crisis) can be cured if 
the Affleck-Dine mechanism for baryogenesis works in the early universe. Furthermore, 
the second one can be solved by raising the mass of the Polonyi field up to O(lOTeV) so 
that the reheating temperature by the (pp decay is larger than O(lMeV). Then, the BBN 
starts after the decay of 0p has completed. 

In order to raise Polonyi mass m^p without raising the SUSY breaking scalar masses in 
the observable sector, several ideas have been proposed. One of the attractive proposal is to 
assume a dynamical SUSY breaking scenario. If the SUSY is broken by some dynamics at 
the energy scale Mj much below the gravitational one, the cut-off scale of the model may 
become 0{Mj). In this case, Kahler potential may contain (non-renormalizable) interaction 
terms which are suppressed by powers of Mf^ instead of M~^. Then, the Polonyi field may 
have a mass much larger than that of the gravitino, and hence decays before the BBN starts. 

Alternative approach is to consider no-scale type supergravity models In some 

class of models with no-scale type Kahler potential, the gravitino mass is not directly related 
to the SUSY breaking masses of the squarks, sleptons and Higgs bosons. In this case, 
the gravitino mass (and the mass of the Polonyi field) can become much larger than the 
electroweak scale without raising the SUSY breaking scalar masses in the observable sector. 

Furthermore, in some class of models with fine tuning, mixing between the Polonyi field 
and the scalar field concerning the fine tuning becomes substantially large In this case, 
the Polonyi field decays through the mixing, and the decay rate can be enhanced without 
increasing the Polonyi mass m^p. If the mixing is sufficiently large, the Polonyi field decays 
before the BBN starts. An example of this type of models is the minimal SUSY SU(5) 
model with a small self-coupling of the 24-Higgs. In this model, the Polonyi field 
decays into the fiavor Higgs doublets through the mixing to 24-Higgs. 

As we have discussed, there exist some approach to the Polonyi problem. However, they 
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still have uncertainties, and it is unclear whether these scenarios really works well. Therefore, 
more efforts are needed in order to check the validity of each scenarios. Finally, it should 
be noted that the cosmological evolution of the Polonyi field may affect the arguments on 
the gravitino problem. In this thesis, however, we assume that the Polonyi field plays no 
significant role in cosmology due to some (unknown) mechanism, and hence we ignore its 
effects in the following chapters. 
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Chapter 6 

Cosmology with unstable gravitino 



As mentioned in the previous chapter, the gravitino may affect cosmology. Especially 
the constraints from the BBN and the present mass density of the universe strongly suggest 
the inflation if there exists the gravitino. In this chapter, we will consider the effects of the 
gravitino on the inflationary universe mainly assuming that the gravitino decays only into 
photon and photino. We will also give the analysis for the case where the gravitino only 
decays into neutrino and sneutrino.Q 



6.1 Gravitino production in the early universe 

We will first calculate the number density of the gravitino after the inflation. Once 
the universe has reheated, gravitinos are reproduced by scattering processes of the thermal 
radiations and they decay with the decay rate of the order of 771^2/ ^pi- Since the interactions 
of gravitino are very weak, the gravitinos cannot be thermalized if the reheating temperature 
Tr is less than 0{Mpi). In this case, the Boltzmann equation for the gravitino number density 
?T,3/2 can be written as 

0^^3/2 , ,^ ,2 "•^3/2 ^3/2 

+ = {^totVrel} U^ad " Je^^' ^ 

where H is the expansion rate of the universe, (■ ■ ■) represents thermal average, rirad the 
number density of the scalar boson in the thermal bath, 

Urad = ^T', (6.2) 

Vrei the relative velocity of the scattering radiations {{vrei) = 1 in our case), and the factor 
1^3/2/ {E3/2) the averaged Lorenz factor. For the radiation dominated universe, the expansion 
^This chapter is based on the work in a coUaboration with M. Kawasaki H, [47||. 
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rate H of the universe is given by 




(6.3) 

where R is the scale factor and A^^, the effective total number of degrees of freedom for 
effectively massless particles. For the particle content of the MSSM, N^:(Tji) ~ 228.75 if Tr 
is much larger than the masses of the superpartners, and N^{T <C IMeV) ~ 3.36. The total 
cross section E^t in the thermal bath is defined as 

^tot = 1^ ^^^y (^(x+y^'iP^+z) , (6.4) 

where cr(^+y^^^+2) is the cross section for the process x + y ^ + z (see Table ^l2|) , and 
rjx = 1 for incoming bosons, rjx = \ for incoming fermions. For the MSSM particle content, 
Stoi is given by 

^tot = ^ {2.50^?(T) + Amgl{T) + 11.78^|(T)} , (6.5) 

where gi, g2 and (73 are the gauge coupling constants of the gauge group U(l)y, SU(2)^ 
and SU(3)(^, respectively. Notice that in high energy scattering processes, the effect of the 
renormalization of the gauge coupling constants should be taken into account. Using the one 
loop /3-function of the MSSM, the solution to the renormalization group equation of gauge 
coupling constants is given by 

,.(r).{,-(,n,)-A_,„(i_)}"'\ (6.6) 

with hi = 11, 62 = 1, ^3 = —3. In this thesis, we use this formula. 

At the time right after the end of the reheating, the first term dominates the right-hand 
side of eq. ( |6.1D since gravitinos have been diluted by the de Sitter expansion of the universe 
during the inflation. As a first step to solve eq.( |6.lD , we assume a "naive" adiabatic expansion 
of the universe; 

RT = const. (6.7) 



Then using the yield variable 13/2 = n3/2/nrad and ignoring the decay contributions, eq.( |6.1| ) 
becomes 

^ = Hf • 

Notice that the right-hand side of this equation is (almost) independent of T, and hence we 
can easily integrate eq.(|6.8|). 
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However, eq.(|6.8D does not give a correct value of F3/2 since the conserved quantity is not 
RT but the entropy per comoving volume; 



R^S = const, 



(6.9) 



with 5* being the entropy density. Therefore, if the number of the gravitino per comoving 
volume is conserved, the yield variables 13/2 for different temperature Ti and T2 are related 



as 



>^3/2(Ti, 



S{Ti)/nrad{Ti) 



NsiT, 



H^3/2(T, 



(6.10) 



where Ns = S/Sq with 5*0 = (27r/45)T^. (The prefactor Ns(Ti)/Ns(T2) is sometimes called 
a dilution factor.) 

Taking this effect into account, the yield of the gravitino at T is given by0 



^ .^x _ NsiT) nrad{TR) {T.totVrel) 
) — — - X 



Ns{Ti 



R) 



H{Tr) 



(6.11) 



For the MSSM particle content, Ns{Tr) ~ 228.75 and Ns{T < IMeV) ~ 3.91. Eq.(|Ol|) 
shows that I3/2 is proportional to Tr, i.e. as the reheating temperature increases the yield 
of the gravitino becomes larger. From eq.( |6.11| ), we can derive a simple fitting formula for 



3/2 1 



Y^I2{T < IMeV) ~ 2.14 x 10" 



-11 



R 



lOiOGeV 



1 -0.0232 In 



R 



lOiOGeV 



(6.12) 



where the logarithmic correction term comes from the renormalization of the gauge coupling 
constants. The difference between the exact formula ( |6.11| ) and the above approximated one 
is within ~ 5% for 10^ GeV < Tr < 10^^ GeV (and within ~ 25% for 10^ GeV < Tr < 10^^ 
GeV). Notice that the gravitino abundance derived here is about four times larger than the 
one obtained in ref. |3^ . 

As the temperature of the universe drops and the Hubble time H^^ approaches 
the decay term becomes the dominant part of the right-hand side of eq.( |6.1| ). Ignoring the 
scattering term, eq. (|6.1|) can be rewritten as 



3/2 



3/2 



dt 



T3/2 



(6.13) 



^In the recent article, Fischlcr |Q have proposed a new mechanism to produce gravitinos in the thermal 
bath. If one adopts his mechanism, the number density of the gravitino becomes much larger than the results 
obtained in this thesis, and hence the constraints obtained below must become more stringent. However, it 
is not clear to us if his new mechanism is relevant. 
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where we have taken m^i2/ (-E'3/2) = 1 since gravitinos are almost at rest. Using eg. ( |6.11| ) as 
a boundary condition, we can solve eq.( |6.13|) and the solution is given by 



"'^^ uraait) Ns{Tn) H{Tn) 



(6.14) 



where the relation between t and T can be obtained by solving eq.( |0|) with eq. (|6.7|) ; 



1 90Af2 „ 

6.2 Radiative decay of the gravitino 

Radiative decays of the gravitino may affect BBN. First, we analyze this effect assuming 
that the gravitino ip^ mainly decays into a photon 7 and a photino 7. 

In order to investigate the photo-dissociation processes, we must know the spectra of the 
high energy photon and electron induced by the gravitino decay. In this section, we will 
derive these spectra by solving the Boltzmann equations numerically. 

Once high energy photons are emitted in the gravitino decay, they induce cascade pro- 
cesses. In order to analyze these processes, we take the following radiative processes into 
account. 

• The high energy photon with energy produces an pair by scattering off the 
background photon if the energy of the background photon is larger than ml/e^. We 
call this process the double photon pair creation. For sufficiently high energy photons, 
this is the dominant process since the cross section or the number density of the 
target is much larger than that of other processes. A numerical calculation shows 
that this process determines the shape of the spectrum of the high energy photon for 

> ml/22T. 

• Below the effective threshold of the double photon pair creation, high energy photons 
lose their energy by the photon-photon scattering. But in the limit of — >^ 0, the total 
cross section for the photon-photon scattering is proportional to and this process 
loses its significance. 

• Finally, photons lose their energy by scattering off charged particles in the thermal 
bath. The dominant processes are pair creation in the nuclei and the Compton scat- 
tering off the thermal electron. 

• Emitted high energy electrons and positrons lose their energy by the inverse Compton 
scattering off the background photon. 

• The source of these cascade processes are the high energy photons emitted in the decay 
of gravitinos. Notice that we only consider the decay channel V'a* — >■ 7 + 7 and hence 
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the energy of the incoming photon e^o is monochromatic. 



The Boltzmann equations for the photon and the electron distribution function and 
fe are given by 



dt 



dt 



+ 



dt 



+ 
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dfe{Ee 



cs 



cs 



dt 



(6.16) 
(6.17) 



IC 



where DP (PP, PC, CS, IC, and DE) represents double photon pair creation (photon-photon 
scattering, pair creation in nuclei, Compton scattering, inverse Compton scattering, and the 
contribution from the gravitino decay). Full details are shown in Appendix p. 

By solving the Boltzmann equations ( |6.16D and (|6.171 ) numerically, we obtain the time 
evolution of the photon and the electron spectra. The typical time evolutions of the photon 



spectrum are shown in Figs. 5.1 



6.3 BBN with high energy photon injection 

BBN is one of great successes of the standard big bang cosmology. It is believed that 
light elements with atomic number less than 7 are produced when the cosmic temperature 
is between IMeV and lOkeV. Theoretical predictions for abundances of light elements are 
excellently in good agreement with the observations if the baryon-to-photon ratio rjB is about 
3 X 10"^° 1|9|. (For a review, see Appendix 0.) 

However, the presence of the gravitino might destroy this success of the BBN. Gravitino 
may have three effects on BBN. First the energy density of the gravitino at T ~ IMeV speeds 
up the cosmic expansion and leads to increase the n/p ratio and hence ^He abundance also 
increases. Second, the radiative decay of gravitino reduces the baryon-to-photon ratio and 
results in too baryon-poor universe. Third, the high energy photons emitted in the decay 
of gravitino destroy the light elements. Among three effects, photo-dissociation by the high 
energy photons is the most dangerous for the gravitino with mass less than ~ ITeV. In the 
following we investigate the photo-dissociation of light elements. 

The high energy photons emitted in the decay of gravitinos lose their energy during mul- 
tiple electro-magnetic processes described in the previous section. Surviving soft photons can 
destroy the light elements (D, T, ^He, ^He) if their energy are greater than the threshold of 
the photo-dissociation reactions. In this thesis, we consider the photo-dissociation reactions 



listed in Table 3J-. For the process D(7,n)p, we use the cross section in analytic form which 
is given in ref . , and the cross sections for other reactions are taken from the experimental 
data. (For references, see Table |6.1|) . We neglect the reactions '^He(7, D)D and ^He(7, 2p 
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(GeV) 

Figure 6.1: Time evolution of the photon spectrum with the case of Tf> = lO^^GeV, and 
(a) 772,3/2 = lOOGeV, and (b) 171^/2 = ITeV. The sohd hues (the dotted hues, the dashed 
hues, and the dotted-dashed hues) correspond to the photon spectra at the time t = O.IT3/2 
(t = T3/2, t = 5t3/2, and t = IOT3/2). 59 



Reaction 



Threshold (MeV) Reference 



D +7 
T +7 ^ 
T +7 -> 
^He +7 
^He +7 
^He +7 
^He +7 
^He +7 



n + p 
n+ D 
p + n + n 

T 

-> ra+^He 
p + n+ D 



2.225 
6.257 
8.482 
5.494 
7.718 
19.815 
20.578 
26.072 



m 

0, n 



m 



Table 6.1: Photo-dissociation reactions. 



2n) since their cross sections are small compared with the other reactions. Furthermore, we 
do not include the photo-dissociation processes for ''Li and ^Be because the data of the cross 
section for ''Be is not available, and hence we cannot predict the abundance of ^Li a part of 
which comes from ^Be. 

The time evolution of the light elements are described by 



"df 
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+ 3HnB 



dt 



dt 



J Ej 

+ E / de 
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J Ej 

-^0(^7) A (^7) 

•/ Ej 

+ Y de^al^^ti^i 



(6.18) 



(6.19) 



(6.20) 



(6.21) 



where ci*^^ is the cross section of the photo-dissociation process i; a + 7 — >• ■ ■ ■, and Ei is 
the threshold energy of reaction i. When the energy of the high energy photon is relatively 
low, i.e. 2MeV ^ e-y ^ 20MeV, D, T and ^He are destroyed and their abundances decrease. 
On the other hand, if the photons have energy high enough to destroy "^He, it seems that 
such high energy photons decrease the abundance of all light elements. However, since D, 
T and ^He are produced by the photo-dissociation of ^He whose abundance is much higher 
than the other elements, their abundances increase or decrease depending on the number 
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density of high energy photon. When the number density of high energy photons with their 
energy greater than ~ 20MeV {i.e. the threshold energy for ^He destruction) is extremely 
high, all light elements are destroyed. But as the photon density becomes lower, there is 
some range of the high energy photon density in which the overproduction of D, T and ^He 
becomes significant. And if the density is sufficiently low, the high energy photon does not 
affect the BBN at all. 

From various observations, the primordial abundances of light elements (D, ^He, ^He) 



are estimated 



as 



0.22 <Yp = 



- 

> 1.8 X 10"^ 



< 0.24, 



nil 



< 1.0 X 10-^ 



(6.22) 

(6.23) 
(6.24) 



where pin^ and pb are the mass densities of ^He and baryon. (For details, see Appendix y.) 
The abundances of light elements modified by the gravitino decay must satisfy the observa- 
tional constraints above. In order to make precise predictions for the abundances of light 
elements, the evolution equations (|6.18|) - ( |6.21|) should be incorporated in the nuclear net- 



work calculation of BBN. Therefore, we modify the Kawano computer code |^ to include 
the photo-dissociation processes. 

From the above arguments it is clear that there are at least three free parameters, i.e. 
the mass of gravitino m3/2, the reheating temperature Tr and the baryon-to-photon ratio 
riB- Furthermore, we also study the case in which the gravitino has other decay channels. 
Here, we do not specify other decay channel. Instead, we introduce another free parameter 
Bry which is the branching ratio for the channel ■j/'^ — 7 + 7. Therefore, we must study 
the effect of gravitino decay on BBN in four dimensional parameter space. However, in the 
following arguments it will be shown that the baryon-to-photon ratio ?7b is not an important 
parameter in the present calculation because the allowed value for r/s is almost the same as 
that in the standard case {i.e. without gravitino). 

First we investigate the photo-dissociation effect when all gravitinos decay into photons 
and photinos {B^ = 1). We take the range of three free parameters as lOGeV < ^3/2 < 
lOTeV, lO^GeV < Tr < lO^^GeV and 10"^° < Vb < 10"^ In this calculation, we assume 
that the photino is massless.0 The contours for the critical abundances of the light elements 
D, (D-|-^He) and ^He in the t]b - Tr plane are shown in Fig. ^]2|- Fig. |6.5| for 7/13/2 = lOGeV, 
lOOGeV, ITeV and lOTeV. For lower reheating temperatures (Tr^ lO^GeV), the number 
density of the gravitino is very low and hence the number density of the induced high energy 
photons is too low to affect the BBN. Therefore, the resultant abundances of light elements 

^Constraints from the photo-dissociation of the hght elements are almost independent of the photino 
mass if the photino is sufHciently lighter than the gravitino. 
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Figure 6.2: Contours for critical abundance of light elements in the ris - plane for 
"^3/2 = lOGeV. The solid line (dashed line, dotted line, dotted-dashed line) represents the 
constraints from overproduction of (D+^He) (overdestruction of D, overdestruction of ^He, 
overproduction of "^He). 
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Figure 6.3: Same as Fig. 3.2 except for = lOOGeV. 
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Figure 6.4: Same as Fig. |6]^ except for m^i2 = ITeV. 
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Figure 6.5: Same as Fig. pl2| except for 7713/2 = 10Te"V. 
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Figure 6.6: Upperbound on Tr as a function of Here, we take = 1. In the region 

above the sohd curve ^He and D are overproduced, the abundance of ^He is less than 0.22 
above the dotted curve and the abundance of D is less than 1.8 x 10~^ above the dashed 
curve. 



are the same as those in the standard BBN. The effect of the photo-dissociation due to 
gravitino decay becomes significant as the reheating temperature increases. 

As seen in Fig. p.2\ - Fig. |6.5| , the allowed range of the baryon-to-photon ratio is almost 
same as that without gravitino for m3/2 ^ ITeV, i.e. a very narrow range around rjB ^ S x 
10~^° is allowed. However, for the cases of ^3/2 ~ ITeV and Tr ~ lO^GeV or 777,3/2 ~ ITeV 
and Tr ~ lO^^GeV , lower values of tib are allowed (Fig. |6.4| ). In this case, the critical photon 
energy (~ m1/22T) for the double photon pair creation process is lower than the threshold 
of photo-dissociation reaction of '^He. Therefore, for Tr ^ lO^^GeV, the abundance of ^He 
is not affected by the gravitino decay. Then, the emitted photons only destroy ^He and D 
whose abundances would be larger than the observational constraints for low baryon density 
if gravitino did not exist. Therefore one sees the narrow allowed band at Tr ~ lO^GeV 
where only a small number of ^He and D are destroyed to satisfy the constraints ( |6.23D and 
{^^. For r^ > IQi^GeV, since a large number of high energy photons are produced even 
above the threshold of double photon pair creation, a part of "^He are destroyed to produce 
■^He and D, which leads to the very narrow allowed region at Tr ~ lO^^GeV. However, 
even in this special case the upper limit of allowed reheating temperature changes very little 
between rjB = 10~^° and t]b ^ 3 x 10~^°. This allows us to fix t]b = 3.0 x 10~^° in deriving 
the upperbound on the reheating temperature. 

The allowed regions that satisfy the observational constraints ( |6.22| ) - ( |6.24| ) are also 
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shown in Fig. EM in the 1713/2 - Tr plane for the case of r^^ = 3 x 10 and Bj = 1. In 



Fig. |6]^ - Fig. ^]5| and Fig. one can see four typical cases depending on Tr and 171^/2. 



• m3/2 ^ ITeV, Tr < WGeV; 

In this case the lifetime of the gravitino is so long that the critical energy for the double 
photon process (~ m1/22T) is higher than the threshold of the photo-dissociation 
reactions for ^He at the decay time of gravitino. Thus ^He is destroyed to produce T, 
^He and D. (Since T becomes ^He by /5-decay, hereafter we denote T and ^He by the 
word "^He".) Since the reheating temperature is not so high, the number density of 
gravitino is not high enough to destroy all the light elements completely. As a result, 
^He and D are produced too much and the abundance of ^He decreases. To avoid 
the overproduction of (^He + D), the reheating temperature should be less than ~ 
(10^ - 109)GeV. 

• m3/2 < ITeV, Tr > lO^^GeV; 

The lifetime is long enough to destroy ^He and the gravitino abundance is so large that 
all the light elements are destroyed since the reheating temperature is high enough. 
This parameter region is strongly excluded by the observation. 

• lTeV<m3/2<3TeV; 

The lifetime becomes shorter as the mass of gravitino increases, and the decay occurs 
when the double photon pair creation process works well. If the cosmic temperature 
at t ~ r3/2 is greater than ~ mg/22ii^4He (where Ei^e ~ 20MeV represents the typical 
threshold energy of ^He destruction processes), ^He abundance is almost unaffected 
by the high energy photons as can be seen in Fig. In this parameter region, the 
overproduction of (D+'^He) cannot occur since ^He is not destroyed. In this case, the 
destruction of D is the most important to set the limit of the reheating temperature. 
This gives the constraint of Tr< {10^ - lO^^^QeV. 

• m3/2^3TeV; 

In this case the gravitinos decay when the temperature of the universe is so high that 
all high energy photons are quickly thermalized by the double photon process before 
they destroy the light elements. Therefore, the effect on the BBN is negligible. Fig. |6.5| 
is an example of this case. The resultant contours for abundances of light elements are 
almost identical as those without the decaying gravitino. 

So far we have assumed that all gravitinos decay into photons and photinos. But if other 
superpartners are also lighter than the gravitino, the decay channels of gravitino increases 
and the branching ratio for the channel ip^ — 7 + 7 becomes less than 1. In this case, 
various decay products affect the evolution of the universe and BBN. In this thesis, instead 
of studying all decay channels, we consider only the 7 + 7 channel with taking the branching 
ratio as another free parameter. With this simplification, the effect of all possible decay 
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products other than photon is not taken into account. Therefore, the resultant constraints 
on the reheating temperature and the mass of gravitino should be taken as the conserva- 
tive constraints since other decay products may destroy more light elements and make the 
constraints more stringent. 

Although we have four free parameters in the present case, the result for = 1 implies 
that the allowed range of Tr and ^3/2 is obtained taking the baryon-to-photon ratio to be 
3 X 10^^'^. Since our main concern is to set the constraints on Tr and m3/2, we can safely 
fix r/B as 3 X 10"^°. 

The constraints for = 0.1 and = 0.01 are shown in Fig. |6.7| and Fig. ^]8| which 
should be compared with Fig. {i.e. the case of -B^ = 1). Since the number density of 
the high energy photons is proportional to i?^, the constraint on the reheating temperature 
becomes less stringent as B^ decreases. In addition, the total lifetime of gravitino is given 
by 

r3/2 = r(V'M^7 + 7)x5^. (6.25) 

Thus the gravitinos decay earlier than for the Bj = 1 case, and hence the constraints from 
(^He + D) overproduction becomes less stringent. This effect can be seen in Fig. |6.7| , where 
the constraint due to the overproduction of ('^He + D) has a cut at 171^/2 — 400GeV compared 
with ~ ITeV for B^ = 1. 

In Fig. |6]^, the contours for the upperbound on reheating temperature are shown in 
the m3/2 - B^ plane. One can see that the stringent constraint on Tr is obtained for 
^^3/2 ~ lOOGeV even if the branching ratio is small. Notice that the actual constraint may 
become more stringent if we take the effects of other decay products into account. 



6.4 BBN with high energy neutrino injection 

As seen in the previous section, the existence of the unstable gravitino which decays 
into a photon and a photino set a stringent upperbound on the reheating temperature Tr. 
However, the constraints might become much weaker when the gravitino decays only into 
weakly interacting particles. In the particle content of the MSSM, the only candidate is the 
decay into a neutrino and a sneutrino. 

Even if the gravitino decays only into a neutrino and a sneutrino, the emitted high en- 
ergy neutrino may scatter off the background neutrino and produce an electron-positron (or 
muon-anti-muon) pair, which then produces many soft photons through electro-magnetic 
cascade processes and destruct light elements. Since the interaction between the high en- 
ergy neutrinos and the background neutrinos is weak, it seems that the destruction of the 
light elements is not efficient. In fact, Gratsias, Scherrer and Spergel showed that the 
constraint is not so stringent for the case where the gravitino decays into a neutrino and 
a sneutrino. However the previous analysis seems to be incomplete in a couple of points. 



First, Gratsias et al. [58| totally neglected the secondary high energy neutrinos which are 



produced by the neutrino-neutrino scattering. The effect of the secondary neutrino may 
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Figure 6.7: Same as Fig. p]0 except for = 0.1. 
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Figure 6.8: Same as Fig. |6]^ except for = 0.01. 
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Figure 6.9: Contours for the upper limits of the reheating temperature in the m3/2 - 
plane. The numbers in the figure denote the limit of the reheating temperature. 



be important for the heavy gravitino case (7723/2 ^ITeV). Second, they only studied the 
case where the destruction of "^He results in the overproduction of (^He + D). However, for 
the heavy gravitino which decays in early stage of the BBN, the destruction of D is more 
important since the electro-magnetic cascade process is so efficient that the energy of soft 
photons becomes less than the threshold of ^He destruction. Furthermore, as pointed out 
in the previous section, the previous estimation of the gravitino production in the reheating 
epoch after the inflation is underestimated. Those effects which are not taken into account 
in ref.|5^ may lead to a more stringent constraint on the reheating temperature. Therefore, 
in this section we reexamine the effects of the decay of the gravitino into a neutrino and a 
sneutrino (7/'^ ^ z/ + z>) with taking all relevant processes into account |]47 |. 

The high energy neutrinos v produced in the gravitino decay scatter off the thermal 
neutrino in the background through the following processes; 



+ yi,h - 


I/j + Z/j, 


+ ^i,h - 




+ ^i,h - 




+ yjf, - 




+ yj,h - 




+ ^ih - 


e" + e+ 



(6.26) 
(6.27) 
(6.28) 
(6.29) 
(6.30) 
(6.31) 
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Figure 6.10: Upperbound on the reheating temperature from BBN and the present mass 
density of the sneutrino. The region above the curves is excluded. 



Z/j + z/, 



i,b 



+ ft 



(6.32) 



where index i and j represent e, /i and r with i ^ j. The primary and secondary high energy 
neutrinos scatter off the background neutrinos and produce charged lepton pairs. Then, 
the charged leptons produced in the processes (|6.31| ) and (|6.32| ) induce the electro-magnetic 
cascade processes. By the same procedure as in the previous section, we obtain the high 
energy photon and the electron spectra. The detailed analysis of the neutrino spectrum is 
given in Appendix p. 

The photo-dissociation of the light elements are analyzed in the same way as in the 
previous section. In the present calculation, there are at least three free parameters, i.e. the 
mass of the gravitino ^3/2, the reheating temperature Tr and the baryon-to-photon ratio 77^. 
However as shown in the previous section, the baryon-to-photon ratio r/s is not important 
parameter because the allowed value for is almost the same as that in the standard case 
[i.e. without gravitino). Therefore, we fix = 3 x 10^^° in the following analysis. 

The allowed regions that satisfy the observational constraints (|6.22|) - (|6.24|) are shown 
in the ^3/2 - plane in Fig. p.lO] . In Fig. |6.1(]| one can see that for the gravitino of 
mass between lOOGeV and ITeV, the overproduction of D and ^He gives the most stringent 
constraint, while the upperbound on the reheating temperature is determined from the 
destruction of D for m3/2 ~ (1 — 3)TeV. Notice that D destruction was not considered in 
the previous work [^. Furthermore, the constraint from (D+'^He) overproduction is more 
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stringent. The reasons why we obtain the more stringent constrain are (i) the gravitino 
abundance is (4 — 5) times larger than the one that the previous authors used, (ii) the 
secondary neutrinos are taken into account in our calculation, and (iii) the gravitino lifetime 
for ip^ ^ u + i) is longer by a factor 2 than the lifetime for -i/^^ — 7 + 7. (In ref.|5^ it is 
presumed that T 3/2(4'^ ^ u + u) = T 3/2(4'^ 7 + 7)0 



6.5 Discussion about hadron injection 

In the previous sections, we have assumed that the gravitino only decays into non- 
hadronic particles, and derived constraints on the reheating temperature. If the gravitino 
produces hadronic decay products, however, the previous constraints may change. In this 
section, we will briefly discuss the effect of such strongly interacting decay products. 

The effects of the hadronic decay products on the BBN are investigated by several au- 
thors ^ and the standard BBN scenario may be affected by the following processes. 

• Injecting high energy particles at the time (1 — 10^)sec cause non-standard p ^ n 
conversion processes. The major effect is to induce p ^ n reactions because there are 
more target protons than target neutrons. The extra neutrons produced in this period 
results in overproduced ^He abundance. 

• Hadronic decay of the gravitino with lifetime T3/2 ^10^ sec causes dissociation processes 
of light elements. Especially even in the case T3/2 ^ lO^sec in which the light element 
photo-dissociation processes are not effective, ^He may be significantly destroyed to 
produce D and ^He. 

In the following, we discuss the above effects. 

The p ^ n conversion induced by hadronic injection was studied by Reno and Seckel |B9 



in details. In the standard BBN scenario, the neutron fraction in the thermal bath changes 
from ~ 0.16 (at T ~ 0.7MeV) to ~ 0.12 (at T ~ O.OSMeV). Hadronic injection at this period 
(O.TMeV ^ T ^ O.OSMeV) affects this evolution. Reno and Seckel claimed that the yield of 
the gravitino Y3/2 should be smaller than 0(10^^^ — 10^^^) in order not to overproduce '^He. 
From this bound, we can estimate that the reheating temperature Tr should be lower than 
O((10^° - 10")GeV) for the case lOTeV < mg/z < lOOTeV. 

If the gravitino lifetime is longer than ~ lOOsec, gravitinos decay after the "deuterium 
bottleneck" breaks. In this case, emitted hadrons induce light element destruction processes, 
which give more stringent constraint than that from ^He overproduction. These processes are 
analyzed in ref.|5^ (for the case T3/2 ^ lO^sec), and in refs. pO|, |6T[] (for the case T3/2 ~ lO'^sec). 

For the case lO^sec^ T3/2 ~ lO^sec, (D-f^He) overproduction induced by "^He destruction 
gives a constraint. According to ref.|]59|, the reheating temperature Tr should be lower than 
O((10^ — 10^^)GeV) if the gravitino mass is given by (l-lO)TeV. 

If the gravitino lifetime is longer than ~ lO^'sec, photo-dissociation of light elements 
becomes effective. In this case, we should take the effects of high energy photons as well as 
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injected hadrons into account. As discussed in refs. |^|, the main effects of the hadrinic 
injection in this case is the destruction of ^He, and the creation of D, ^He, ^Li and ^Li. 
According to refs. [^0[ the hadronic branching ratio must be very small in order not to 
overproduce (D+'^He), and hence the constraint is almost the same as in the case = 1, 
except for the case T3/2 ~ 10^~^sec. In a small parameter region ~ 10^~^sec, one may 
avoid the severe constraint even in the case 5^ = 1 by assuming the non-vanishing hadronic 
branching ratio, since the effect of photo-dissociation of D is compensated by the supply 
of D through the ^He destruction processes. However in this case, ^Li and ''Li seem to be 



overproduced. Furthermore, as discussed in ref.|^0[ some uncertainties exist in the analysis 
of refs. pU| |BT| . The most important source of uncertaintie is in the experimental data of 
hadron scattering processes, especially those concerning Li. Therefore, it is unclear to us if 
the scenario proposed in refs. EO, BI| works well. 



6.6 Other constraints 

In the previous sections, we have considered the constraints from the photo-dissociation 
of light elements. But as we have seen, if the mass of the gravitino is larger than a few TeV, 
the gravitino decay does not induce light element photo-dissociation and no constraints has 
been obtained. In this case, we must consider other effects of the gravitino. 

If we consider the present mass density of the LSP produced by the gravitino decay, we 
can get the upperbound on the reheating temperature. In SUSY models with i?-invariance 
(which is an usual assumption), the LSP (in the previous cases, photino or sneutrino) is 
stable. Thus the LSPs produced by the decay of gravitinos survive until today, and they 
contribute to the energy density of the present universe. Since one gravitino produces one 
LSP, we can get the present number density of the LSP ulsp as 

riLSP = Yy2{T « IMeV) x ^T^^ow, (6-33) 

where T^ow — 2.7K is the present temperature of the universe.0 The density parameter of 
the LSP 

i^LSP = , (6.34) 

pc 

can be easily calculated, where ttilsp is the mass of the LSP, pc — 8.1 x 10~'^'^h'^GeV'^ is the 



In calculating eq.( 6.33 ), we have ignored the effect of the pair annihilation of the LSP (7 + 7 —> / + /, or 
D+D* — * f + f )- For the pair annihilation processes, the cross section is roughly estimated as crv ^ oiw /''^lsp 
(where aw is the coupling factor and v is the relative velocity), and the condition for sufficiently large pair 
annihilation rate {ulsp'^v ^ H, with ulsp being the number density of the LSP) reduces to riLSp/n-y ~ 10~^ 



(with Yj^sp being the yield variable for the LSP), which is less stringent than the constraint (5.35). However, 
if the LSP mass is small enough to hit the pole of the Z°-boson propagator, the constraint (3.35) may become 



weaker, and the constraint from the "^He overproduction may become significant for the large gravitino mass 
case (m3/2 ^ 3TeV). 
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Figure 6.11: Upperbound of the reheating temperature. Dashed hue represents the constraint 
from the present mass density of photino. Dotted-dashed hne represents the upperbound 
requiring ^He < 0.24. Constraints from D photo-dissociation is also shown by dotted hne. 



critical density of the universe and h is the Hubble parameter in units of lOOkm/sec/Mpc. If 
we impose ^Ilsp < 1 in order not to overdose the universe, the upperbound on the reheating 
temperature is given by 



Tr < 2.7 X 10 



11 / rriLSP \ ^ 
UOOGeV/ 



GeV, 



(6.35) 



where we have ignored the logarithmic correction term of T^tot- 

To set the upperbound on the reheating temperature, we need to know the mass of the 
LSP. First let us consider the case where the LSP is photino. In this case, if one assumes the 
gaugino-mass unification condition, the lower limit of the mass of photino is 18.4GeV ||62|| . 
Then we can get the following upperbound on the reheating temperature; 



Tr < 1.5 X lO^'rGeV (if the LSP is photino). 



(6.36) 



This constraint is shown in Fig. |6.11| . Notice that this bound is independent of the gravitino 
mass and branching ratio. If the LSP is sneutrino, the present limit on sneutrino mass 



41.8GeV 1 24] sets the upperbound on the reheating temperature; 



Tr < 6.6 X GeV (if the LSP is sneutrino) 



(6.37) 
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which is also shown in Fig. |6.10. 



Another important constraint comes from the effect on the cosmic expansion rate at the 
BBN. As mentioned before, if the density of gravitino in the nucleosynthesis epoch becomes 
higher, the expansion of the universe increases, which leads to more abundance of "^He. We 
study this effect by using the modified Kawano code and show the result in Fig. |6.11| . In 
the calculation, we take tib = 2.8 x 10^^^ and r„ = 887 sec (where r„ = (889 ± 2.1)sec is 
the neutron lifetime [^) so that the the predicted ^He abundance is minimized without 
conflicting the observational constraints on other light elements. The resultant upperbound 
on the reheating temperature is approximately given by 



Tr<2x lO'^GeV ( ) , (6.38) 



for m3/2 > lTeV.0 This bound is also shown in Fig. |6.11| by the solid line. Notice that this 
bound is derived irrespective of the species of the LSP. 



^For the case < ITeV, primordial ^He is destroyed by the high energy photons induced by the decay 
of the gravitino, and hence the constraint ( 6.38| ) becomes insignificant. However, if the gravitino decays into 
a neutrino and a sneutrino, this is not the case. 
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Chapter 7 

Cosmology with stable gravitino 



In the previous chapter, we have studied cosmology with unstable gravitino and derived 
constraints on the reheating temperature after inflation. However, there is another possibility 
that the gravitino is stable and hence it is the LSP. In this case, the constraints obtained in 
the previous chapter are not appropriate. In this chapter, we study cosmological constraints 
when the gravitino is the LSP.Q 



7.1 Constraints from the mass density of the universe 

In the case of the stable gravitino, the Boltzmann equation for the gravitino number 
density 71,3/2 is given by 

^ + 3Hny, = + E -I^r, (7.1) 

where nj is the number density of the superparticle i, Tj the decay rate of i into its super- 
partner and a gravitino, and m--/ (E:-) the Lorentz factor. Notice that the first term in the 
right-hand side of eq. ( [7. 1|) represents the contribution from the scattering processes off ther- 
mal radiations, while the second one is that from the decay of superparticles into a gravitino 
and some ordinary particles.0 

Here we comment on TilK'^^ . If the gravitino mass is small compared with the typical mass 
splitting in the matter sector, interactions of the helicity ±| modes of the gravitino become 
stronger than that of the helicity ±| modes. From this fact, we conclude that the hehcity 
±1 gravitino production processes are the most significant in this case. As we can see in 
eq.( |4.73| ), the interactions of the gravitino to gauge multiplets are described by dimension- 



five operators and those to chiral multiplets by dimension-four ones. Therefore, the former 



^This chapter is based on the work in a collaboration with H. Murayama and M. Yamaguchi [ p^ . 

^In the case of heavy unstable gravitino analyzed in the previous chapter, we have ignored the contribution 
from the decay of the superparticles (or its inverse processes). In fact, these contributions are significant only 
for the case of extremely light gravitino (7713/2 ^ 10~^GeV, as we will see later), and hence our approximation 
is justified. 
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dominates the gravitino production at high temperatures. For the MSSM particle content, 
Siof^ is given by 



247rm2/2M2 



{2.Uglml, + 9.16glml, + 26mglml:,) , (7.2) 



where a{x + y ^ ip + ^) is the hehcity ±^ gravitino production cross section (see Table] 
^G3 ^ ™Gi are the gauge fermion masses, and g^ - gi are the gauge coupling constants. 

Using the yield variable Yx = nx/nrad and the condition of "naive" adiabatic expansion 
(|0|), eg. 1|) is rewritten as 

dT HT ^ HT ' ^ ' ' 

i 

Integrating this equation from Tr to T (Tr ^ T) and taking the effects of the dilution factor 
{Ns{T) / Ns(Tr)) into account, we get 

Yy,{T) = X {YscattiT) + F,eea,(T)} . (7.4) 

with 

i scatt — rr/rp \ ) I'-^j 

^ decay T ^ {E') H 

From eq.( [r.5|) , we can see that Y scatt is proportional to the reheating temperature Tr. On 
the other hand, Ydecay is almost independent of Tr as far as the reheating temperature is 
higher than the masses of the superparticles. Notice that effects of the gravitino annihilation 
processes cannot be ignored in eq . ( [7. 1|) if the Yg^att in eq. (|7.5|) or the Ydecay in eq . ([71^ ) 
becomes 0(1). We take these effects into account by the following simple method; in the 
case where the sum of the right hand sides of eq. (|7.5| ) and eq. (|7.6|) , which is the yield of 
the thermal helicity zb^ gravitino, is larger than |, we consider that the gravitinos are 
thermalized and take Y scatt + Ydecay = §• 

Using eq.( [7.5|) and eq.( [7.6|) , we estimate the present gravitino number density and de- 
termine the upperbound on the reheating temperature by using the closure limit. In our 
numerical calculation, we take all the squark and slepton masses to be 1 TeV and the GUT 
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Figure 7.1: Cosmological constraints on the gravitino mass ms/2 and the reheating temper- 
ature Tr in the framework of the MSSM when the gravitino is the LSP. We take all the 
squark and slepton masses to be 1 TeV, rriGi = rnj^sp = 50GeV and the GUT relation on 
the gauge fermion masses is assumed. The solid line denotes the upperbound on the reheat- 
ing temperature from the closure limit. The dotted line is the upperbound on the gravitino 
mass from light element photo-dissociation. We assume the NSP relic density as given in 
eq.( |7.11| ), and require the NSP lifetime to be shorter than 2.6 x 10^ sec. 



relation on the gaugino masses; 

3 mci 



(7.7) 



5 gf 92 gi 

Energy density of the gravitino P3/2 is given by 

P3/2(T) = m3/2F3/2(T)n^ad(T), (7.8) 

and we determine the upperbound on Tr by the condition 



P3/2{Tnow) < Pc, (7.9) 

where Tnqw - 2.7K and Pc - 8.1 x IQ-^'^h^ GeV^ (0.4 < /i < 1) is the critical density of the 
present universe. 

The result is shown in Fig. |7.1| . Here we take m^i = rriNSP = 50 GeV with rriNSP being 
the mass of the next-to-the-lightest superparticle (NSP), and to get a conservative bound we 
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take h = 1. For 171^/2 ^ 10 GeV, the upperbound on Tr is approximately proportional to 
7713/2. When 2 x lO'^ GeV<m3/2< 10"^ GeV, the upperbound on Tr is around 0(100 GeV). 
And for the very light gravitino case (7/13/2 ~ 2 x 10~^ GeV), there is no constraint on Tr. 
These features can be understood in the following way. When the gravitino mass is larger 
than about 10"'* GeV, interactions are so weak that decaying processes cannot produce 
sufficient gravitinos to overdose the universe. Therefore, it is the scattering process that is 
important to estimate the number density of the gravitino. In this case. 



^ . NsjTNow) V90C(3)M ^ 

Y3/2{TnOw) = ^^^y^^ i^totVrel) , (7.10) 



from eg. ( [7. 51) . Combining eq.( [7.10|) with eq.( [7.9|) , we get the upperbound on the reheating 



temperature, which is approximately proportional to the gravitino mass. On the other hand, 
if 2 X 10~^ GeV^7n3/2^10"'^ GeV, the decay processes become significant. In this case, P3/2 
is larger than pc unless the reheating temperature is smaller than the squark and slepton 
masses. Therefore, it is necessary to lower the reheating temperature below the squark and 
slepton mass scale in order not to overdose the universe. And when 7/13/2 ^ 2 x 10^^ GeV, 
the gravitino mass is so small that P3/2 cannot exceed pc even if the gravitino is thermalized. 

7.2 Constraint from BBN 

Next, let us consider the constraint from the light element photo-dissociation. If a decay 
of a heavy particle produces high energy photons after the primordial nucleosynthesis, we 
must require that these photons do not change the abundance of the light elements. Here we 
consider the decay of the NSP. Since the gravitino is the LSP, the NSP can decay only into 
a gravitino and some ordinary particles through the supergravity interaction. Therefore, the 
NSP have much longer lifetime than other superparticles and may affect the predictions of 
the BBN. 

If the NSP were stable, it would survive until today. Its relic density in this case has 



been calculated []63| , |64| , [65[] . For the neutralinos, in a wide range of parameters, the relic 
density is larger than 10^^ to the critical one. This relic density can be translated into 
i^nspYnsp > 5.0 X 10^^^ GeV,0 where rriNSp and Y^sp are the mass and yield of the NSP. 
In the following analysis, we conservatively take 

mNSpYNSP = 5.0 X 10"^^ GeV, (7.11) 

and assume that the NSP decay produces high energy photons. Photo-dissociation of light 
elements by the NSP decay can be analyzed in the same way as in the case of unstable 
gravitino. Detailed analysis shows that the energy density of the NSP as large as the one 
given in eq. ( [7.11| ) will overproduce (^He + D) unless the lifetime of the NSP is shorter than 
•^It is plausible that this bound is also valid when a slepton or a chargino is the lightest. 
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about 2.6 X lO^sec. (See Fig. - Fig. |8.3| in the next chapter.) Therefore, we impose 



1 ^NSP 



ASnmhoM'^ 1 KniNSpJ 



"3/2^ 



< 2.6 X lO^'sec. 



(7.12) 



Here, we assume that the NSP is the U(l)y gauge fermion (bino) and used eq.( [4.74| ) for the 
decay rate of the NSP.0 The right hand side of eq. ( [7.12| ) strongly depends on the NSP mass 
and especially when the NSP mass is small, a severe upperbound on the gravitino mass is 
obtained. The bound we obtained is m3/2 < 1-2 GeV (6.7 GeV, 244.4 GeV, 711.0 GeV) for 
rriNSP = 50 GeV (100 GeV, 500 GeV, 1 TeV). This bound for the case of rriNSP = 50 GeV 
is also shown in Fig. [7.1| . If the reheating temperature is sufficiently low compared to the 
NSP mass, the NSP is not produced significantly and the above constraint can be avoided. 
Furthermore, the constraint is not significant when the NSP decay does not produce sufficient 
photons. This is the case if the NSP is a sneutrino. 



7.3 Remarks 

Before closing this chapter, we would like to consider phenomenological implications of 
the gravitino LSP. As we discussed previously, the lightest among the superpartners of the 
standard model particles is no longer stable. Therefore, this particle (NSP) can be charged 
or even colored. In a supergravity model with a no-scale like Kahler potential squarks 
and sleptons are massless at the tree-level. Although they acquire their masses from gaugino 
masses through radiative corrections, a right-handed charged slepton tends to be the lightest. 
Indeed with the GUT relation of the gaugino masses, the mass of the right-handed slepton is 
very close to the U(l)y gaugino mass. Detailed analysis showed that a neutral superparticle 
is the lightest only when its mass is lighter than ~ 150 GeV |6^. The gravitino LSP 



can cure this difficulty of the fascinating class of supergravity model which might be derived 
from superstring compactification It can also solve the conflict |]TD|, |7T], |72| between 

the constraints from the mass density of the neutralino LSPs and those from the proton 
decay in the minimal SUSY SU(5) GUT, because the former gets meaningless when there 
is a superparticle lighter than the neutralino. In accelerator experiments, the most drastic 
change occurs when the NSP is charged. Then, we have a chance to flnd a spectacular track 
of the charged NSP in a detector. 

Finally, we we comment on the case where the gravitino is extremely light (7/13/2 -C 
1 keV) and interacts strongly. In this case, the gravitino decouples from the thermal bath 
at temperature well below the electroweak scale. The standard nucleosynthesis scenario 
constrains the number of the species of the neutrino-like light particles to be smaller than 



^If the bino is the NSP, it decays into a gravitino and a photon, and into a gravitino and a Z'^-boson. But 
when the bino is hghter than the Z°-boson, the latter decay channel is forbidden kinematically and the decay 



rate of the bino is sin^ 9w — 0.23 times smaller than the value of eq.(4.74). For the case rriMSP — 50GeV, 
we have considered this effect. 
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3.3 1^.0 Therefore, in order to get a sufficient dilution, the gravitino must decouple before 
T ~ 200 MeV. After the QCD phase transition, I + 1 ^ ip + ip scatterings become the most 
important in the thermalization processes of the gravitino, where / denotes e, /i, Ue, or 
Vri and Tits anti-particle. Comparing the gravitino production rate of these processes with 
the expansion rate of the universe, we get 

1-13 n^^T ( '"^i 



™3/.>10-"GeVx^^^J. (7.13) 

where is the slepton mass.Q Remember that the interactions of the light gravitino are 
proportional to 'm^^2 as we have seen in Chapter ^, and hence the lowerbound on the gravitino 
mass is obtained here. 



^In a recent paper, Kernan and Krauss [T^ has claimed that the number of the species of neutrinos N^, 
should be smaller than 3.04. In this thesis, however, we use N^, < 3.3 in order to derive a conservative 
constraint. 

^One can also consider the constraints from the collider experiments ^ or from the cooling of the 
red supergiant star Jt^. These arguments give similar bound on the gravitino mass. 
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Chapter 8 

Conclusions and discussion 



8.1 Summary of conclusions 

We summarize the effects of the massive gravitino on the inflationary universe. If the 
gravitino is unstable and decays only into a photon and a photino, stringent constraints 
on the reheating temperature are obtained. In the case where the gravitino mass m^i2 
is less than ~ ITeV, (D +^ He) overproduction due to the photo-dissociation of ^He set 
the upperbound on the reheating temperature Tr, Tr ^ 10^~^GeV. Constraints in the case 
ITeV ^ m3/2 ^ 3TeV, the photo-dissociation of D gives the upperbound on the reheating 
temperature, Tr ^ 10^~^^GeV. If the gravitino mass is heavier than ~ 3TeV, the reheating 
temperature is constrained to be lower than ~ lO^^GeV not to overdose the universe. These 
constraints are shown in Fig. 
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(8.1) 
(8.2) 
(8.3) 
(8.4) 

If most of the decay products of the gravitino are charged particles, above constraints 
are expected to be valid since the spectrum of the high energy photon is fixed mainly by the 
total amount of energy injection. An exception is the case where the gravitino only decays 
into a neutrino and a sneutrino. In this case, the constraints are less stringent than in the 
case of high energy photon injection, and the constraints are given by 

Tr < lO^^GeV (m3/2<100GeV), (8.5) 
Tr < lO^^-^^GeV (100GeV<m3/2<5TeV), (8.6) 
Tr < lO^^QgV (5TeV<m3/2<10TeV), (8.7) 



which are also shown in Fig. |6.10 
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If the gravitino is light and stable, interactions of helicity ±| gravitino become strong 
and a more stringent bound on the reheating temperature is derived from the closure limit; 



Notice that this constraint is not appropriate for the case ^3/2 ^ IkeV. Furthermore, the 
effects of the NSP decay on BBN set the upperbound on the gravitino mass. The detailed 
value of the upperbound depends on the mass of the NSP. For example, the gravitino mass 
larger than ~ IGeV is excluded for the case of tjinsp = 50GeV. 



8.2 Discussion 



We have derived constraints on the reheating temperature after inflation by assuming 
the existence of the massive gravitino field. As we have seen, cosmological arguments give 
stringent upperbound on the reheating temperature. 

Some comments on our results are in order. When the gravitino is unstable and its 
mass is smaller than ~ ITeV, the upperbound on the reheating temperature is given by 
(10^ — 10^)GeV. This constraint seems to be very stringent since such a low reheating 
temperature requires very small decay rate of the inflaton field. For example, in a chaotic 
inflation with a inflaton whose interactions are suppressed by M~^, the decay rate of 
the inflaton is expected to be Finf ~ mf^^/M^i with minf being the inflaton mass, and hence 
the reheating temperature is estimated as 



Tn ~ O.ly/FinfAfp; ~ lO^GeV, (8.9) 

requiring that the inflaton field should produce the density perturbations observed by COBE 
(minf ~ lO^^GeV) |7^. Notice that if the interaction of teh inflaton becomes stronger, the 
decay rate becomes larger and hence a higher reheating temperature is obtained. Thus, in 
the case that the gravitino is lighter than ~ ITeV (which is favored from the naturalness 
point of view), we have to adopt an inflaton with extremely small decay rate. 

Let us compare our results with those in other literatures. The gravitino number density 
and the photon spectrum are the essential quantity for the analysis. Our number density of 
gravitino produced in the reheating epochs of the inflationary universe is about four times 
larger than one given in ref. [3^. Since the authors of ref.[0 neglected the interactions 



between the gravitino and chiral multiplets (which is the second term in eq. ([4.58|) ), they 



might underestimate the total cross section for the production of gravitino. All previous 
works concerning the gravitino problem were based on the gravitino number density given 



in ref. [36 . Therefore our constraints are more stringent than others. 



Furthermore, our photon spectrum is different from that in ref.|40| as shown in Fig,B.2 



The spectrum adopted by ref. [EH] has more power to destroy light elements above threshold 
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Figure 8.1: Contours for critical abundance of light elements in the tx - ejoBjYxo plane for 
eo = lOTeV. The solid line (dashed line, dotted line, dotted-dashed line) is a bound from 
overproduction of (D+^He) (overdestruction of D, overdestruction of '^He, and overproduc- 
tion of ^He). In deriving the mass density of X, we take = 1 and mx = 2e^o- 



for the photon-photon scattering and less power below the threshold. In refs. [|75|, [Tq] , the 
Compton scattering process was not taken into account in calculating the photon spectrum 
which the authors of ref. used to derive a fitting formula for the high energy photon 
spectrum. Therefore, it is expected that the difference comes mainly from the neglect of 
Compton scattering off the thermal electron, which is the most dominant process for the 
relatively low energy photons. 

Finally, we discuss applications of the cosmological arguments to other cases. If some 
exotic particle decays radiatively with a lifetime longer than ~ Isec, the standard BBN 
scenario may be affected. The constraints we have used in the case of massive gravitino can 
be applied to other cases of exotic particle X. Using the similar method as in the gravitino 
case, we derive constraints from the BBN on radiatively decaying weakly interacting particles. 
The results are given in Fig. p.l| - Fig. ^.3| , where tx is the lifetime of X, Yxo the yield of 
X before X decays, e^o the energy of the primary photon, and 5^ the branching ratio of 
radiative decay channel. The evolution of the yield variable Yx = nx/urad of X is given 
by Yx{t) = lxoe~*/^^. Notice that we assume that the emitted photon is monochromatic, 
and that the number of photons emitted from one decay process is one. As we can see, the 
constraint on the combination e^oB^Yxo is almost the same as in the cases of e^o = lOGeV 
and e^o = lOTeV since the spectrum of the energetic photon depends almost only on the 
total amount of energy injection, while for the case of e-^o = lOMeV the constraint from 



82 



10' 



14 r I I I I I I I I I 



10' 10^ 10" 10^ 10^ 10^ 10^ 10^ 10^° 10^^ 10^' 



Xx (sec) 



Figure 8.2: Same as Fig. except for e^o = lOGeV. 
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Figure 8.3: Same as Fig. 



|0 



except for e^o = 10 MeV. 



(D +^ He) is negligible since the energy of primary photon is below the threshold of ^He 
destruction. 

Constraints from the present mass density of the universe may give another information 
on exotic particles. Especially in the SUSY models, the present mass density of the relic 
LSP is calculated in detail, and constraints on SUSY parameters are derived. Furthermore, 
particles with lifetime longer than ~ 10^°sec may distort the cosmic microwave background. 

These arguments have been applied to the LSP in SUSY models, the Polonyi field which 
is responsible for SUSY breaking the massive neutrino the axion |]S2|, ^ 



and so on as well as the gravitino. Thus, cosmological considerations provide us with a great 
insight into exotic particles which we cannot investigate by any collider experiments because 
of the weakness of their interactions. 
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Appendix A 
Notations 



A.l Conventions 

We use the metric (in flat space-time) 

(7^, = diag(l, -1,-1,-1), (A.l) 
and for totally anti-symmetric tensor (in flat space-time), we use the convention €0123 = — 1- 

A.2 Two component notation 



For two component spinor, we essentially follow the notation used in ref. [26], except for 
the convention of the metric, 
cr-matrices are defined as 

^0 = ^0= (J ?)' ^^ = -^i=(? J)' 

a2 = -a2=(° ^^ = -^^=(j -l)- (^-2) 

For the anti-symmetric symbols e^/j and e°^, we use the following conventions; 

= -e^^ = = -ei2 = 1, (A.3) 
ell ^ ^22 ^ ^ ^ Q ^^ 4) 

Notice that e^^ and e"^ are related as eap^^^ = ^2- By using e°^, cx^ and are related as 
From these cr-matrices, the generators of the Lorentz group is the spinor representation 
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are obtained 

CTiiu ^ ^{o-^ai, - ai,aij) , = -{a^a^ - a^a^) . (A.6) 

A. 3 Four component notation 

7-matrices obey the following commutation relation; 

[7m>7J=2^^,. (A.7) 
From 7° - 7"^, 75 matrix is defined as 

75 = ^7"7'7'7'- (A.8) 
7-matrices are transposed by a charge conjugation matrix C; 

C-'^^C = -7;. (A.9) 
Charge conjugation matrix C satisfies the following identities; 

C^C=1, (A.IO) 

= -C. (A.ll) 

By using C, charge conjugation of four component spinor is defined as 

V'^ = Ci^^. (A.12) 

Chiral representation 

Four component spinor (in chiral representation) is constructed from two component 
spinor ^q, and 77^ as 

V^-(<,ed)- (A.13) 
In this representation, 7-matrices are given by 

and charge conjugation matrix C is given by 

C = i7Y = C*. (A.15) 
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Dirac representation 

7-matrices in Dirac representation are related to those in chiral representation through 
unitary transformation; 

^^chiral) _ [/^Pi™c)f;t^ (A. 16) 

with 

Exphcit form of the 7-matrices in Dirac representation are given by 

/^(J_°0. ^-(-./o')' '^-^^^ 

and charge conjugation matrix C is obtained as 

C = i7Y = <^*- (A.19) 
In Dirac representation, solution to the Dirac equation (in momentum space); 

(i.-m)M(p,s) = 0, (A.20) 
with p'^ = (po, |p| sin 9 cos 0, |p| sin ^ sin 0, |p| cos 9) takes the following form; 



where p^p*^ — w?, n = p/|p|, and 

r.=-,i) ^ / cos e/2 \ ^ / -e^W^ sin ^^/2 \ 

^ - 1^ e^'^/2sin^/2 j' ^ e'^/2^os^/2 J" 

The above spinor ii(p, s) is a eigenstate of hehcity nS; 

(nS)M(p, s) = s u(p, s), (A.23) 

where S is given by 
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Appendix B 
Photon spectrum 



In order to investigate photo-dissociation processes, we have to calculate a photon spec- 
trum induced by decay of gravitino. In this appendix, we will write down Boltzmann equa- 
tions which determine the high energy photon spectrum with the case of high energy photon 
injection, and that of high energy neutrino injection. We will also solve them numerically 
and show the shape of the spectrum. 



B.l Boltzmann equations 

In calculating photon spectrum (with high energy photon injection), we take the following 
processes into account; 

• double photon pair creation, 

• photon-photon scattering, 

• pair creation in nuclei, 

• Compton scattering off thermal electron, 

• inverse Compton scattering off background photon, 

• radiative decay of the gravitinos. 

The Boltzmann equations for these cascade processes are formally given by 
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+ 



dt 



+ 



(9/7 (e^ 



pp 



IC 



+ 
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+ 



DE 



dt 
dt 



+ 
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(B.2) 
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where terms with the index DP (PP, PC, CS, IC, DE, and EXP) represents the contribution 
from the double photon pair creation process (photon-photon scattering, pair creation in 
nuclei, Compton scattering, inverse Compton scattering, contribution from the gravitino 
decay, and the effects of the expansion of the universe). Notice that in the main part of this 
thesis, the expansion terms and the decay term in the Boltzmann equation for the electron 
distribution function are ignored. Below, we will see contributions from each processes in 
detail. 



Double photon pair creation [ 7 + 7 ^ + e ] 

For high energy photon whose energy is larger than ~ m?^/22T, double photon pair 
creation is the most dominant process. 

The total cross section ctdp for the double photon pair creation process is given by 



^DP m = \^rl (1 - P^) I (3 - P') In 1±| - 2/? (2 - P^) 

where rg is the classical radius of electron which is given by 

a 



(B.3) 



:b.4) 



with a ^ 1/137 being the fine structure constant, and (3 is the electron (or positron) velocity 
in the center of mass frame. Using this formula, one can write down (9/^/9t)|DP as 



a/^(e. 



dt 



o ,2 ^7 1^7 J 



DP 



7 '-'y 



de^ — fy{e^) / ds saopiP) 



Ami 



:b.5) 



/3=^l-(4m2/s) 



The spectrum of final state electron and positron is obtained in ref. [^, and {dfe/dt)\-DP 
is given by 



dt 



-irr^m^ / ae. 



DP 



7 S 
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A (^7) ^/ TP p \ 
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where function G{Ee, e^, e^) is given by 

v\2 



G{Ee, e^, e^) 



EeE: 
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:b.7) 



with E'^ = e^ + e^ — Ee, and L represents the distribution function of the background photon 
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at temperature T, 



f.,(c,) = 4 X . (B.8) 

Phot on- phot on scattering [7 + 7^7 + 7] 

If the photon energy is below the effective threshold of the double photon pair creation, 
photon-photon scattering process becomes significant. This process is analyzed in ref.|l76 
and for e'^ ^ 0{m\/T), {df^/dt)\pp is given by 

35584 o o _« /-"^ , „ . . o L /e^\ 



dt 



pp 101257r 
1946 



df (e' ) 35584 f°° f e' /"^ 

/ / - 



506257r' 

For a larger value of e'^, we cannot use this formula. But for high energy photons, photon- 
photon scattering is not significant because double photon pair creation determines the shape 
of the photon spectrum. Therefore, instead of using the exact formula, we take ml/T as a 
cutoff scale of {df^/dt)\p-p, i.e. for < ml/T we use eq. (p.9|) and for e'^ > m^/T we take 



df.ie' > ml/T) 



= 0. (B.IO) 

pp 



dt 

Notice that we have checked the cutoff dependence of spectra is negligible. 

Pair creation in nuclei [7 + A^^e+ + e +A/^] 

Scattering off the electric field around nucleon, high energy photon can produce electron 
positron pair if the photon energy is larger than 2me. Denoting total cross section of this 
process cxpc, {df^/dt)\-^p is given by 



•9/7 (£7) 



dt 



= -nArcrpc(e^)/^(e^), (B.ll) 

NP 



where tin is the nucleon number density. For cxpc, we use the approximate formula derived 
by Maximon [^]. For near the threshold (e-^ < 4me), the approximate formula is given 
by 



o-pc(e7)U<4 = —Z'^arj 



3 '\kj 
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where 



(B.13) 



and Z is the charge of nuclei. For large (e^ > 4me), <Jpc are expanded in the parameter 
k~^, which is given by 
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Differential cross section for this process dapc/dE^ is given in ref . \fiT\\ ; 
dupc 
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where 
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In 
In 



E^E_ + J9+J9_ + m; 



E^E_ — p+p^ + m 
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(B.16) 
(B.17) 

(B.18) 



with E^ {E+) being the energy of electron (positron) in final state. By using this formula, 
idfe/dt)\t^P is given by 



dfe{Ee 



dt 



NP 



rfe.^/,(e,). 



' dE, 



(B.19) 
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Compton scattering [ 7 + e ^ 7 + e ] 

Compton scattering is one of the processes by which high energy photons lose their 
energy. Since the photo-dissociation of hght elements occurs when the temperature drops 
below ~ O.lMeV, we can consider the thermal electrons to be almost at rest. Using the total 
and differential cross sections at the electron rest frame acs and dacs/dEe, one can derive 



dt 
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cs 



n, 
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(B.20) 
(B.21) 



where n^. is the number density of background electron. Explicit forms of acs and dacs / dE^ 
are given as 



o'cs 

dacs{e'^,e^) 
de' 
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(B.23) 



with 



X = 



s — ml 2e^ 



(B.24) 



Inverse Compton scattering [ + 7 ^ + 7 ] 

Formula for the inverse Compton process is given by Jones [Q, and {df /dt)\ic is given 

by 



5/7 (ei 



dt 
dfe{E: 
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where function -F(e^, i^e, e^) is given by 

Ffe.,E.,e-.)L . = 2glng + (l + 2g)(l-g) + — (B.27) 



-7, -^ej '-7/l0>(j>l 
F(e-y, i?e, £7) I otherwise 



2 (1 - T,q) 



0, 



(B.28) 
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with 



Radiative decay of the gravitino [ ^ 7 + 7 ] 

Source of the non-thermal photon and electron spectra is radiative decay of gravitino. 
Since gravitinos are almost at rest when they decay and we only consider two body decay 
process, incoming high energy photons have monochromatic energy e^o? which is given by 



£70 



m2/2 - m? 



2m: 



■3/2 



(B.29) 



In this case, (9/^/9t)|DE can be written as 

•9/7 (£7 



dt 



B^5 {ey 



DE 



^3/2 
T3/2 ' 



(B.30) 



where is the branching ratio for the process V'^t ~^ 7 + 7- If the gravitino decays into an 
electron and a selectron, (i9/e/i9t)|DE cannot be ignored. Denoting the branching ratio for 
this decay process Be, the formula for (9/e/9t)|DE is given by 



dfe{Ee 



dt 



BJ (Ee-E, 
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3/2 



DE 



T3/2 



(B.31) 



where Eeo is the energy of the injecting electron. 



B.2 Spectrum with high energy photon injection 

Having obtained the explicit forms for the Boltzmann equations, we solve them in this 
section. In deriving high energy photon and electron spectra, we have to use numerical 
method since the Boltzmann equations and ( p.2|) are very much complicated. Before 

solving them numerically, we explain our approach to the Boltzmann equations p9 |. 

At first, we classify the right-hand side of eq. (|B.l|) and eq. (|B.2| ) into the "outgoing" parts 
and "incoming" ones. Below, we ignore the expansion terms (9/-y/(9t)|EXP and {dfe/dt)\-Exp 
in the Boltzmann equations ( [B.l| ) and (p3.2D since the expansion rate of the universe is much 



smaller than the scattering rates of electro- magnetic processes. Then, eq. and eq. (p.2| ) 
can be written as 

^■^"^'"^ = -r,(e,;T)/,(e,) + /,,iN(6,), (B.32) 



dt 

dfe{E, 

dt 



-Te{Ee- T)fe{Ee) + /e,IN(^e), (B.33) 
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where "incoming" terms /^jn and /e,iN can be obtained as the sums of the contributions 
from the decay terms and the functionals of the distribution functions of the photon and the 
electron; 



/7,IN(e7) 



fe,m{Ee) 



POO 



dt 

+ dEiK,^,{e,,Ei;T)ME'^), (B.34) 

dfe{Ee) 



DE 



+ / de'^Ke,^{Ee,e'^;T)f^{e^ 

J 



dt 

roc 

+ dE:K,AE,,E:-T)UE:). (B.SS) 

The exphcit forms of and Kab -B = 7, e) can be obtained from the full details of the 
Boltzmann equations given in the previous section. 

As mentioned before, the strategy we use here is to calculate a stationary solution to the 
Boltzmann equations ( [B.l| ) and ( |B.2| ), which obey 



•^T. = ' = 0. B.36 
dt dt ^ ^ 

In our approach, the distribution functions f-y and fe of the photon and the electron can be 
formally given as 

f.i^.) = ^?^y (B.37) 
fe{E.) = (B.38) 

The important thing is that the "incoming" terms f-y^ni^-y) and fe,w{,Ef.) depend only 
on /^(e') and /e(e') with e' > e^,Ee. Therefore, if the distribution functions /^(e^) and 
fe{Ef,) with e^,Ee > e (and the source terms) are known, we can obtain f.y{e) and /e(e) 
from eq.( p.37|) and eq. (p.38| ). By using this fact, we solve the Boltzmann equations with 
monochromatic high energy photon injection. As one will see, extensions to the cases with 
non-monochromatic sources or the high energy electron source are trivial. 

In our numerical calculations, we use the mesh points Ei {0 < i < A^mosh + 1) with 
Eq = Cmin and Ejy^^^^^ = eo, where emin is some minimum energy we are concerning (which 
we take emin = IMeV) and eo is the energy of the primary injecting photons. The energy 
range between emin and eo is divided logarithmically, i. e. the i-th mesh point E^ is given by 

^i = emin^ . (B.39) 

\ ^min / 

By using the fact that the incoming photons into the mesh point E^^^^i^ are supplied only 
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by the photon source, incoming terms for i = A^mesh is given by 



A,in(^^_J ^ n^x-^, (B.40) 

* mcsn 



feMEN^.sJ ^ 0, (B.41) 

where is the production rate of monochromatic photons, and Aj = (-Ej+i — i?j_i)/2. Com- 
bining eq. (|B.40| ) and eg . (16.411) with eq. (|B.37|) and eq.( |B.38|) , we can obtain the distribution 



functions at the A^mesh-th mesh point /^{Ejy^^^^^J and /e(-EAr^g^^). 

Next we determine the distribution functions for lower energy region. Essentially, "in- 
coming" terms for the i-th mesh point are derived from the distribution functions f-yiEj) 
and fe{Ej) with j > i. Discretizing eq.( [B.34|) and eq. (|B.35| ), "incoming" terms are (approx- 



imately) given by 

f,ME^) - Y.^AK,,{E„E,)f,{Ej)+K,,{E„E,)UE,)}, (B.42) 

j>i 

feMEi) ^ Y.^AKe^{E,,E,)f^{Ej)+K,,{Ei,E,)UE,)}, (B.43) 



from which we can obtain f.y{Ei) and fe{Ei) by using eq.( p.37|) and eq. (p.38| ) .p] In a way 



explained above, we calculate the photon and electron distribution functions f-y{Ei) and 
fe{Ei) at each mash points from i = A^mesh — 1 to z = in order. 

For each T and e^O) "we calculate the spectra /-y(e^) and fe{Ee) by solving eq.( p.l| ) and 
eq. (|B.2|) numerically. Typical spectra are shown in Fig. in which we show the case with 



e-yo = lOOGeV and lOTeV, T = lOOkeV, IkeV, lOeV, and the incoming flux of the high energy 
photon is normalized to be 



dt 



5(e^ - e^o) GeV^ (B.44) 



DE 



Notice that eq. (|B.l|) and eq. (p.2|) are linear in and fe- Therefore, once the solution /^,ref 



is obtained with some reference value of the decay term (<9/y/5t)|DE,ref, we can reconstruct 
the photon spectrum for arbitrary value of (5/^/9t)|DE with temperature T and the incident 
photon energy e^o fixed; 

M^l) = Ureii^j) X fr., ,^,^| • (B.45) 

[OJ-y/ Crj|DE,ref 

The behaviors of the photon spectrum can be understood in the following way. In the 
region ml/22T, the photon number density is extremely suppressed since the rate of 
double photon pair creation process is very large. Just below this threshold value, the shape 



iln fact, eq.( |B.40| ) - eq. (|B.43D receive corrections of the order of A, which is due to the discretization of 



the energy range. In our numerical calculations, we take these corrections into account. 
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Figure B.l: Typical spectra of photon (the sohd hnes) and electron (the dotted lines). We 
take the temperature of the background photon to be T = lOOkeV, IkeV, lOeV, and the 
energy of the incoming high energy photon is (a) lOOGeV and (b) lOTeV. Normalization 

nf fVip initial nVinfnn is crivpn 1= ^ \' ( f) f ( /= ^/ f)i\\ -r^-r, = S( t: — t: ^\ fipV^ 



6^0 = 10 TeV 



Temperature 


Plow iVi^^GeV^ Ppp iVppGeV^ 


1 eV 
10 eV 
100 eV 


-1.57 1.6 X 10* -5.10 6.9 x 10"^* 
-1.34 5.4 X 10* -5.20 6.0 x 10"^* 
-1.22 1.7x10^ -4.84 1.1 x 10^^^ 


e-yo = 1 TeV 


Temperature 


Plow iVio^GeV^ Ppp iVppGeV^ 


1 eV 
10 eV 
100 eV 


-1.56 1.4 X 10* -5.07 6.2 x lO'^* 
-1.34 4.9 X 10* -5.17 5.5 x lO'^* 
-1.22 1.4 X 10^ -4.79 1.0 x lO^^^ 


e^o = 100 GeV 


Temperature 


Plow iViowGeV^ Ppp iVppGeV^ 


1 eV 
10 eV 
100 eV 


-1.56 1.4x10* -5.01 5.7x10-1* 
-1.33 4.7 x 10* -5.15 5.3 x lO'^* 
-1.22 1.3 X 10^ -4.74 1.1 x lO^^^ 


e^o = 10 GeV 


Temperature 


Plow iViowGeV^ Ppp iVppGeV^ 


1 eV 
10 eV 
100 eV 


-1.33 4.5 X 10* -5.12 5.5 x 10"^* 
-1.22 1.3 X 10^ -4.77 9.6 x lO^^* 



Table B.l: Plow, -^low, -Ppp and A^pp for the cases of T = 1 eV, 10 eV, 100 eV, and e^o = 
10 TeV, ITeV, lOOGeV, lOGeV. Here, we take iV^esh = 401. 



of the photon spectrum is determined by the photon-photon scattering process, and if the 
photon energy is sufficiently small, the Compton scattering with the thermal electron is the 
dominant process for photons. 

As one can see in Figs. |B.1|, both the photon spectrum for sufficiently low energy region 



and that for the photon-photon scattering region obey power-law spectrum; f-y{e-y) oc e:^. 
We fit the photon spectrum at the sufficiently low energy region as 

Pin \ f T / 
VG 

and that for the photon-photon scattering region as 
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with Pin being the total amount of the energy injection from the gravitino decay; 



Pin 



dt 



(B.48) 



DE 



Notice that the amphtude of the spectrum is proportional to the mean free time of the 
photon. In low energy region the mean free time is determined by Compton scattering and 
depends on the background temperature as ~ T~^. For the photon-photon scattering region, 
the amplitude is proportional to ~ T^^ since the cross section depends on ~ T^. For the 
cases T = leV, lOeV, lOOeV and eo = lOTeV, ITeV, lOOGeV, lOGeV, we calculate the 
fitting parameters Piow, -Ppp, ^low and A^pp, and the results are shown in Table |B.l| .p| These 
variables slightly depend on the background temperature, while their dependence on the 
initial photon energy e^o is insignificant. 

In Fig. B.2 , we compare our photon spectrum with the results of the simple fitting formula 
used in ref . , which is given byQ 



where 
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^3/2 rfnij(e^ 



T3/2 



de.. 



24^2 



0<e^<emax/2 



55 ^ 



^-3/2 



de^ 



55 



max^7 1 



duE^e^ 



de^ 



(B.49) 



(B.50) 
(B.51) 
(B.52) 



with 



22T' 



(B.53) 



As one can see, not only the absolute value but also the form of the spectrum differs 
between them. The fitting formula in ref.]^ is derived from the numerical results given in 
refs.[^, in which, however, the effect of the Compton scattering is not taken into account. 
Our results indicate that the number of Compton scattering events is comparable to that 
of the inverse Compton events for such low energy region, since the number density of the 
high energy electron is extremely smaller than that of high energy photon. Therefore, the 



^Notice that the formulae ( p3.46D and ( |B.47 ) are relevant for the cases when the initial photon energy e^o 
is much larger than the effective threshold of the double photon pair creation process. 

^Although the photon spectrum is not explicitly given in ref. pof , we obtain it from their photon pro- 
duction spectrum divided by the cross section for Compton scattering. Since the cross section has energy 
dependence, the resultant spectrum (cx e^*^'^) becomes softer than that for photon production (c>c e^^'^)- 
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Figure B.2: Photon spectrum derived from the fitting formula used in Elhs et al. |^0| (solid 
line) is compared with our result with = lOOGeV (dashed hne). We take the temperature 
of the background photon to be lOOeV and the normalization of the incoming flux is the 
same as Fig. [B.l. 



deformation of the photon spectrum by Compton scattering is expected below the threshold 
of the photon-photon scattering. 

Before closing section, we should comment on the case with high energy electron sources. 
By numerical calculations, we have checked that the spectra in the case with high energy 
electron injection are almost the same as those with photon injection if the background 
temperature and the total amount of the energy injection are fixed. Numerically, their 
differences are at most 0(10%). Therefore, the formulae given in eq.( [B.46| ) and eg. ( |B.47] ) 
with Table |B.1| are well approximated ones if one redifines the piN as 



Pin 



df- 



dt 



DE 



dt 



(B.54) 



DE 



B.3 Spectrum with high energy neutrino injection 

Next we will consider the case of high energy neutrino injection. In this case, the emitted 
high energy neutrinos may scatter off the background neutrinos and produce an e^e~ (or 
p+p") pairs, which then produces many soft photons through electro- magnetic interactions. 
Since the rates of neutrino-neutrino scattering processes are not large enough, we cannot 
ignore the effect of the expansion of the universe. Therefore we calculate high energy photon 
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^i + Vi 
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16 sin^ 6' VI/ 
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Table B.2: Coefficients a - d for each processes. Index i and j (with i ^ j) represent the 
generation, i/j is the background neutrino of i-th generation, if is the charged lepton of i-th 
generation (in our case, or fi"^), and 9^ is the Weinberg angle. 



spectrum in two steps; ffist we determine the time evolution of the distribution function of 
high energy neutrinos, and then we calculate the photon (and the electron) spectrum by 
regarding the high energy neutrinos as sources of high energy e~^e~ (and ^Ji~) pairs. The 
procedures for the second step are essentially equal to those given in the previous section, 
and in this section, we will explain the first step. 

Let us begin by deriving the Boltzmann equations which determines the high energy 
neutrino spectra. Injecting high energy neutrinos scatter off the thermal neutrinos in the 
following processes; 

i^i + i^i,b + i^i, (B.55) 

+ ^i,b —> + ^i, (B.56) 
i^i + i^i,b ^j + ^j: (B.57) 
+ ^j,b + ^j, (B.58) 

i^i + i^j,b — > i^i + i^j, (B.59) 
z/j + z/j_fe + e^, (B.60) 

i^i + i^i,b (B.61) 

where index i and j are the generation indices with i ^ j, and represents and yU^. All 
the amplitude squared I^Mp in these reactions take the form given by 

\M\^ = 32Gl {a{pp'Y + b{pqY + c{pq'f + dm^pp')} , (B.62) 

where Gp = {l/^v^) ~ 1.17 x lO'^GeV"^ is the Fermi constant (with v ~ 246GeV), p 
and p' the initial momenta of high energy neutrino and background neutrino, q and q' the 
final momenta, m the mass of the fermion in final state, and the coefficients a - d depend 
on individual reaction. Coefficients for each processes are given in Table |B.2| . 

First, we will consider the neutrino scattering processes ( |B.55D - ( |B.59D . Let us derive 
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a formula for the energy distribution of scattered neutrino for the case where a high energy 
neutrino with energy runs through a region filled with neutrino gas with temperature T^. 
For this purpose, it is convenient to consider the momenta p- q' in the center-of-mass frame. 
In the center-of-mass frame, we parametrize p - g' by using parameters C, t], ( {0 < ^,r] < tt, 
and < C < 27r) as 



p^^ = 00(1, 0, sine, cosO, (B.63) 



p'^^ = 00(1, 0, -sine, -cosO, (B.64) 
q^!^ — uj{l, sinrycosC, sin77sine, cost]), (B.65) 
9cm = '^(1) — sinr^cosC, — sin?7sinC, —cost]). (B.66) 



Then the differential cross section is given by 

da 



^—Gl {aipp'f + h{pqf + cipc^f + dm\pp')} , (B.67) 



d cos rjdC 27i'^s 
with 

s = {p + p'f = Auj^. (B.68) 

Boosting this system to z-direction with velocity we obtain 

p'' = a;(7(l + /3cosO, 0, sin^, 7(008^ + , (B.69) 
p'^ = a;(7(l-/3cosO, 0, sin^, 7(- cos^ + , (B.70) 

with 7~^ = VI — P"^- We identify these momenta as those in the comoving frame. Then the 
energy of initial high energy neutrino Ei, and that of thermal neutrino E^ are given by 

E, = a;7(l + /?cos0, (B.71) 
E^ = u;j{l- pcos^). (B.72) 

Furthermore, the angle 6 between p and p' (in the comoving frame) can be obtained from 
the following relation; 

2 - 23"^ 

cos^^l- ^ (B.73) 
1 — cos e 

and the energy of the neutrino in final state E'^ is given by 

El = q° ^ u-f (1 + (3 cos ri). (B.74) 
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Using these identities and integrating with differential cross section ( [B.67| ) becomes 
da 



dE' 



327r 



-1 



8a + (3 — cos ^ — cos i] — A cos ^ cos r/ + 3 cos ^ cos 1 
+c {?> — cos^^ — cos^?7 + 4cos^ cos?7 + 3cos^^ cos^r^^ |. (B.75) 

Now we are ready to calculate the number density of scattered neutrino with energy E'^ 
- [El + dEl) in unit time. For this purpose, let us consider the process in which incident 
high energy neutrino with energy E^, - [E^, + dE^j) scatters off the thermal neutrino with 
energy E^ - {E^ + dE^) and relative angle 9 - {9 + d9). For this process, number density of 
the target (thermal) neutrino is given by 



(target) = - f,y{E^)dE^ sin 9d9. 
where fu{Eu) is the distribution function of the background neutrino; 



27r2 e-E^/T'^ + r 



(B.76) 



with Ti, being the neutrino temperature. Furthermore, relative velocity of the two neutrino 
is (1 — cos^^) and hence one can obtain the incident flux as 

(flux) = {1- cos 9) UE,)dE,. 

Then one can obtain the contribution to the time derivative of the neutrino distribution 
function as 



dfAK) 



dE' 



(target) x (flux) x 



da 
'dE[, 



1287r 



X 



1 



1 - ^ - 

Gip I dE„ I dE,j „ - fy{Ey)fy{Ey 



EyE^ 



ds (E- 



I (8a + 36 + 3c) — 4(6 — c) cos ^ cos rj 

+ (6 + c) {?) COS^ ^ COS^ 7] — COS^ ^ — COS^ Tj^ |, 



(B.77) 



with 



i?in = Ey + Ey. 



(B.78) 



Notice that the angles ^ and r] can be represented by the variables in the comoving frame 
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through the following relations; 



cos^ = 
cosr^ = 

Define 

Hp (^Ey, E^j 
then eq. (pT^ ) can be rewritten as 



E,-EJ [El-s 



{2El - Ein) (eI - s 



-1/2 
1/2 



(B.79) 
(B.80) 



ds (^El 



-p/2 



(B.81) 



dE' 



X {(8a + 36 + 3c)Hi (e^, E^) 
-(6 + c) {E,-E,y H^{e,,E,) 
-{b + c) {2El-E,^fH^(E,,E,) 
-4(6 - c) (e, - E,) {2EI - E,^) H, [e,, E, 
+3(6 + c) (e, - K)\2EI - E,^f H, [e,, E,) }. 

After some complicated calculations, the above equation becomes very simple form; 



(B.82) 



dUK) 



dEl 



4 r - E - - 

—Gl / dE, / dE,-^UE,)UE,) 
37r JEi J Ef, 



x\aEl + h{E,- Elf + cE^ 



+ ^aE,E, + h {E, - El) E, - ^cElE, 
+^{a + 6b + c)EA. 



(B.83) 



In the following analysis, we use the approximation E,,El ^ E,. With this approximation, 
the above equation becomes 



dUK, 



4 

3^ 

X 



Gl / dE,— iaEl + h{E,- El) + cE^H UE,) 
J El Ef, I- J 

dE,Ej,{E,). (B.84) 



Next we will calculate how many neutrinos scatter off the background neutrino in unit 
time. As is the similar method used above, contribution to the time derivative of the neutrino 
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distribution function from this effect can be obtained as 



BE,. 



4 / 1 1 \ r°° 
-—Gli^a + -b + -cjEME,)J^ dE^EME,), (B.85) 



wliere cr(s) is the total cross section obtained from the amphtude ( B.62 ). Notice that the 
condition for the neutrino number conservation is reahzed; 



dUE, 



OEjy 



oEi, 



(B.86) 



unless the effects of inerastic channels (i/ + z/ — e"*" + e~, /i^ + are taken into account. 

Effects of the charged lepton pair creation process can be taken into account in the same 
way, and the contribution to the time derivative of the neutrino distribution function is given 
by 



dUE, 



dEl 



with 



4 \ ^ Am 



X 



2 ^ 1/2 



(B.87) 



-Anr In 



E^E^ (SElEl - 2m^E^E^ - 3m^ 
2 {4 - (4mVE,^,) }'^' E,E, + AE.E, - 



2m? 



(B.88) 



h 



2 U - 



4m' 



2 ^ 1/2 



EyEy 



-2m In 



2 {4 - (4mVE,E,) }'^' + 4E,^, - 2m' 
2m^ 



(B.89) 



Coefficients for the charged lepton production processes are given in Table [B.2 . 

By using the above formulae, one can obtain the Boltzmann equation describing the 
evolution of the spectra for the high energy neutrinos; 



dt 



I "pff"^ ^^"^ XI {0'in,ijEl + bin,ij {Ey - E'^Y + Cin;ijE'Y^ 
J E' Hi,, 
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X / dE.EME, 
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Stt V 3 3 / JO 



+ 



dt 



+ 



hi^i— >e-+e+ 



9t 



1 



n^/25{El - m3/2/2) 



6r3/2 



(B.90) 



where if is the expansion rate of the universe, and the coefficients a - c are given by 



Q-out — 4, bout — 0, Cout — 13, 

(^in,ii 6, bifida 9, Cifi^ii 11, 



a 



in,t] 



1) binij 1, 



(B.91) 
(B.92) 
(B.93) 



The formula for charged lepton pair creation are given in eq.( [B.87|) . 

We solve the Boltzmann equations for the distribution function of z/g, i/^ and Ur numer- 
ically. In Figs [B.3| - |B.5| , we shown the time evolution of the distribution function of the 
electron neutrino for the cases of 7713/2 = lOOGeV, ITeV, and lOTeV. We have also checked 
that the difference among three types of neutrinos are very small since the differences only 
comes from the charged lepton pair creation processes which are sub-dominant compared 
with the neutrino- neutrino scatterings. In Fig. B.6 , the spectra of three types of neutrinos 
at T = leV are shown for 777,3/2 = lOOGeV, ITeV, and lOTeV. (The time evolutions of the 
photon spectra are also shown in Fig. [B.3| - Fig. [B.5| .) 

By regarding these high energy neutrinos as sources of high energy charged leptons {i.e. 
e^e~ and fi~^fi~), we calculate the high energy photon spectrum induced by the gravitino 
decay into a neutrino and a sneutrino. The Boltzmann equations for the photon and the 
electron distribution function and /e are given by 
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+ 



dt 



+ 



•9/7 (£7) 



pp 



dt 



+ 
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dfe{E, 



cs 



(B.94) 



cs 



dt 



IC 



NEU 



dt 



(B.95) 



EXP 



where NEU represents the contribution from the u - u scatterings. In our numerical calcu- 
lations, we neglect the effects of the expansion of the universe as in the previous case. In 
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Figure B.3: Time evolution of the distribution function of (a) the electron neutrino, and 
(b) the photon for the case m3/2 = lOOGeV. The solid curve (dotted curve, dashed curve, 
and dotted-dashed curve) represents the spedlQTim at the time t — 10^^X3/2 (T3/2, IOT3/2, and 
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Figure B.5: Same as Fig. p.3| except for m3/2 = lOTeV. 
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Figure B.6: Spectra of high energy electron-neutrino (sohd curve), muon-neutrino (dotted 
curve) and tau-neutrino (dashed curve) at T = leV for 7713/2 = lOOGeV, ITeV and lOTeV. 
The gravitino abundance is normahzed as 13/2 = e~*'''^3/2_ 



solving eg. ( p.94|) and eq. (|B.95|) , we assume 



fe{E) 



dt 



-E 



NEU 



dt 



+ 



*e-+e+ 



dt 



(B.96) 



which respects the energy conservation although the actual energy distribution may be dif- 
ferent. This assumption is adequate for our purpose since the photon (electron) spectrum 
is determined almost only by the total amount of the energy injection . Furthermore we 
treat muons (and anti-muons) as electrons with the same energy and neglect the contribution 
from tau leptons whose creation rate is much smaller than the other charged leptons. Full 
details for other terms are shown in the previous section. 
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Appendix C 

Big-bang nucleosynthesis 



Along with the existence of the cosmic microwave background, big-bang nucleosynthesis 
(BBN) is one of the most important predictions of the big-bang cosmology. As we will see 
later, observations of the abundance of light nuclei with atomic number less than 7 are in 
good agreements with theoretical predictions. In this appendix, we will briefly review the 
theoretical framework of BBN, and compare its results to the observations. 

C.l Theoretical framework 

BBN occurs when the temperature of the universe drops below ~ IMeV. In this section, 
we will briefly review the BBN scenario by following the thermal history of the universe of 
temperature about IMeV. (More details, see ref. |]75|.) 

T > IMeV 

When the temperature of the universe is much larger than IMeV, nuclear statistical 
equilibrium (NSE) among the light nuclei is established. Especially, proton (p) and neutron 
(n) are converted to each other in the following weak interaction processes; 

n <-> p + e~ + I7e, 
n + Uf, ^ p + e~ , 
n + p + Fg. 

Since the rate of p n conversion is sufficiently large, chemical potential of e, z/g, p and n 
are related as 

/Un + /i:. = /ip + /ie, (C.l) 

where fix represents the chemical potential of particle X. 

Here, we comment on the magnitude of the chemical potentials of leptons. Since the 
universe is charge neutral, charge density of electron is equal to the number density of 
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proton if T is much smaller than the muon mass. If the temperature is high enough, electron 
can be regarded as a relativistic particle. In this case, charge density is 0{fj,eT^) and we 
obtain 



T 



^!!£ = !^^o(io-^°). (C.2) 



T > IMeV rij 



In order to estimate the neutrino chemical potential, we must know the lepton number 
density of the universe. In the early universe when the temperature is higher than the 
electroweak scale, baryon number and lepton number are converted to each other through 
the sphaleron transition. If this is true, lepton number of the universe is as the same order 
as the baryon number. From this fact, we adopt fJ^u/T ~ O(10~^°). 

In kinetic equilibrium, number density of nuclei with atomic number A and charge Z is 
given by 

UA = 9A J exp j , (C.3) 

where qa represents the internal degrees of freedom. Especially number density of p and n 
can be written as 

If the nucleon with the atomic number A and charge Z can be made out of Z protons and 
{A — Z) neutrons rapidly enough, ha is related to /ip and /i„ in the following way; 

llA = Zfip + {A- Z)fin. (C.6) 

Using eq.( |C.3|) - eq. ( |C.6| ), ua can be expressed as 

. 27r \3{^-i)/2 

UA = gAA'1^2'^ [£^) nlnt' exp (E^T) , (C.7) 

where Ba is binding energy which is defined as 

Ba = ZrUp + {A — Z)mn — mA, (C.8) 

and for simplicity all the "nucleon masses" (m^, m„ and m^/A) in the pref actor are replaced 
by the common mass rriN since their differences are not important. Numerical values of Ba 
and qa for some nuclei are shown in Table |(J.1| . 
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Ba{McY) 


9A 


D 


2.22 


3 




6.92 


2 


^He 


7.72 


2 


^He 


28.3 


1 




92.2 


1 



Table C.l: The binding energy and internal degrees of freedom of light nuclei. 



For a later convenience, we define the mass fraction of species A; 



Xa = (C.9) 



with Ub being the baryon number density; 



Ub ^ np-\-nn-\-^AnA- (C.IO) 

A,Z 



Then we can get 



Xa = yA{C(3)^-V(l-^)/22(3^-5)/2|^5/2(3./^^)3(^-l)/2 

xvi-'X^Xt' exp {Ba/T) , (C.ll) 

where rjB = riB/n^f is the baryon-to-photon ratio. Notice that r]B is related to the present 
baryonic density parameter as 

r]B = — ^ 2.68 X lO-^Qfi/i^. (C.12) 

Since is proportional to r)B~^, mass fraction of species with large atomic number is 
extremely small if they are in chemical equilibrium. For example, with T — lOMeV, 

Xp ~ 0.5, 
Xn ^ 0.5, 

^ 10-'\ 
•v-23 



X3He ^ 10-^^ 



~ 10-3^ 



O.SMeV 5; T ^ IMeV 



Once the temperature of the universe drops below ~lMeV, situation changes. When 
the temperature becomes Tp — 0.8MeV, p^n conversion rate becomes smaller than the 
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expansion rate of the universe. At this freeze out temperature Tp, ratio of Up and n„ is given 
by 



Hp 



T=Tf 



exp ( j ^ -. (C.13) 



Below Tpi number density of neutron decreases mainly through the free neutron decay 
process. But the decay rate of neutron is smaller than the expansion rate and hence, mass 
fraction of neutron is almost unchanged. On the other hand, nucleon production rate is still 
sufficiently large, number density of light nuclei (like D, ^He, ^He, ■ ■ ■) takes the NSE value. 
At this stage mass fractions of light nuclei are as follows; 

Xp ~ 1/7, 
X„ ~ 6/7 



10^ 



12 
23 



XsHe ^ 10 
X4He ^ 10-2^ 



T < O.SMeV 



When the temperature becomes ~0.3MeV, NSE value of X4jje approaches to 1. This 
fact suggests that ^He can be synthesized effectively at r<0.3MeV. But for T>0.1MeV, 
D, which is a source of ^He, is easily destroyed through the scattering processes off the 
background photon. 

Once the temperature drops below ~0.1MeV, D's are effectively produced and they are 
rapidly translated into ^He through the following processes; 

D(D,n)^He(D,p)^He, 

D(D,p)3H(D,n)^He, 

D(D,7)^He. 

Through these processes, almost all of neutrons in the universe are trapped in ^He and 
primordial abundance of D and ^He becomes much less than that of ^He. Order of the 
primordial mass fraction of ^He, Yp, can be estimated to be 



_ 4n4He 

J-p — 

riB 



2n„ 



riB 



~ 0.25, (C.14) 

T~0.1MeV 



where we have used X„|r^o.iMeV ~0.125. (Notice that Yp given in eq.( |C.14| ) is overestimated 
and the result of numerical calculation shows that Yp ~0.24.) 

Finally, we will comment on the nucleosynthesis beyond '^He. Since there is no tightly 
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bounded isotope with A=5 and 8, nucleus which is heavier than ^He is scarcely synthesized. 
One important exception is ^Li which can be synthesized in the following two processes; 



^Hem7)^Li, 
^He(3He,7)^Be(e-,z/e)^Li. 



C.2 Numerical results 

In order to calculate the primordial abundance of light elements, we use numerical 
method. The following results are obtained by using the FORTRAN code written by 



Kawano 57 



Before looking at the numerical results, let us consider about the basic parameters of 
BEN. Essentially, primordial abundances of light elements depend only on one cosmological 
parameter, rjB, however predicted abundances of light nuclei are also affected by uncertain- 
ties come from the uncertainties of the fundamental theory (especially, by those of weak 
interaction) . 

Neutron lifetime r„ is the significant parameter of BEN. In order to calculate the primor- 
dial abundance of ^He, we must accurately determine neutron freeze out temperature Tp, 
which depends on weak interaction rates. These rates are calculated from only one matrix 
element which is also related to neutron lifetime r„, and hence experimental data of neutron 
lifetime is very important for calculating the primordial abundances. We will see that the 
predicted abundance of ^He decreases as neutron lifetime r„ increases. This can be under- 
stood in the following way. A larger value of r„ implies smaller weak interaction rates, and 
this leads to a lower freeze out temperature Tp. Therefore, neutron lifetime increases, ratio 
of Hn to Hp at Tp decreases and hence a smaller value of ^He abundance is expected. 

In Fig. |C.1|, we show the results of the numerical calculations of EEN. For a neutron 



lifetime r„, we use r„ = 891.2sec and r„ = 887.0sec. As mentioned before, for a larger value 
of r„, a smaller value of the primordial ^He abundance is obtained. Primordial abundances 
also depend on the baryon-to-photon ratio rjB- As one can see, predicted abundance of ^He 
is a increasing function of t]b while those of D and ^He decreases as tib increases, and ''Li 
abundance takes its minimum value at ~ 3 x 10"^*'. 

In Fig. |(J.1| , we also plot the predicted abundance of ^He in a model with four neutrino 
species. As one can see, if we add extra neutrino species, predicted abundance of "^He 
increases. This can be understood in the following way. If the extra neutrino exists, the 
expansion rate of the universe at T ~ IMeV gets higher due to the energy density of extra 
neutrino, and hence the neutron freeze out temperature Tp gets higher. In general, if there 
exists some exotic particle whose energy density at T ~ IMeV is not negligible, it speeds up 
the cosmic expansion rate and raise the neutron freeze out temperature. In this case more 
^He is synthesized since lager amount of relic neutron is expected, this argument constrains 
the properties of exotic particles. For example, existence if additional light neutrino species 
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Figure C.l: Predicted abundances of (a) ^He, (b) D and ^He, and (c) ^Li. For a neutron 
lifetime t„, we use r„ = 891.2sec (upper line) and — 887.0sec (lower line). Solid lines 
corresponds to the results with three neutrino species, and for ''He, predicted abundance 
with four neutrino species are also shown in dashed lines. Dotted lines represent bounds on 
riB obtained from observations. 
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conflicts with observational data of ^He abundance. 



C.3 Observations 



In order to see the validity of the theoretical predictions of big-bang nucleosynthesis, we 
must know the actual values of the primordial abundances of light elements. But it is not 
a simple task because the abundances of light elements evolve with time. Especially some 
nuclei are produced or destroyed in stars, which changes the abundances of light elements. 
Therefore we must reconstruct primordial abundances from the observational data. In this 
section, we will derive the constraints on the primordial abundances of D, '^He, "^He and ^Li, 
and compare them with the theoretical predictions. The arguments and observational data 
used in this section mainly follows ref . . 



D and ^He 

We start with discussing the primordial abundance of D. Since the binding energy of D 
is very small (~2.23MeV), D is easily destroyed by (7,p) reaction at 6.0 x lO^K, and it 
is hard to produce D after BBN. Therefore most of D's in the present (or pre-solar) universe 
are expected to be produced by the big-bang nucleosynthesis. In the following we regard the 
observed pre-solar D abundance as a lowerbound on the primordial abundance of D. 

Let us reconstruct pre-solar D abundance from the observational data. Pre-solar value 
of (D+^He) abundance can be found in gas-rich meteorites; 



3.84 X 10"* < 



D + =^He 

me 



< 4.22 X 10" 



(C.15) 



preQ 



On the other hand, carbonaceous chondrites, which is one of the most primitive solar system 
materials, provide us a pre-solar ^He abundance; 



1.48 X 10"^ < 



^He 



< 1.56 X 10"^ 



(C.16) 



preQ 



In order to normalize the abundance of the light elements by hydrogen abundance, we use 
the standard solar model prediction on ^He abundance; 



0.09 < 



and we obtain 



1.8 X 10"^ < 



^He 
IT 



< 0.11, 



(C.17) 



preQ 



< 3.3 X 10" 



preQ 



(C.18) 
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1.3 X 10"^ < 
3.3 X 10"^ < 



^He 
IT 



< 1.8 X 10"^ 



pre© 

D + ^He 



H 



< 4.9 X 10"^ 



(C.19) 
(C.20) 



preQ 



As mentioned before, we expect the primordial abundance of D is larger than the pre-solar 
value since the abundance of D decreases with the galactic evolution. Therefore, pre-solar 
value of D abundance requires that 







D- 




> 




1. 


p 


.1 



> 1.8 X 10" 



pre© 



(C.21) 



Next we will consider the upperbound. Since some of D's are destroyed, it is hard 
to estimate the primordial D abundance without uncertainty. But D is burned to ^He, 
abundance of the sum (D+^He) is less uncertain. Therefore we sill see the upperbound on 
this quantity below. 

Once D is trapped in stars, it either survives or burns to ^He. Furthermore some '^He in 
stars are destroyed. Therefore using the survival fraction of ^He in stars, gsue, primordial 
abundance of (D+^He) can be written as 



D + ^He 
H 



pre© 



D + ^He 
H 



pre© 



'He 

J preQ 
+ {9^1 - 1 



'He 
IT 



(C.22) 



- preQ 



Combining eq. (|C.22| ) with eq.( |C.19|) and eq. (|C.20|) , and taking (yfSHe > 0.25, we can get 

"D + 3He 



H 



< 1 X 10" 



(C.23) 



^He 



In our universe, "^He is the most abundant element next to hydrogen, but some of 
them have non-primordial origin, i.e. some ^He's originate to stars. In order to see the 
uncontaminated primordial abundance of ^He, we shall use the observational data taken from 
the old environments which are not affected by galactic evolution. Here, the primordial ^He 
abundance has been estimated from the observation in metal-poor extragalactic Hnregion. 



In Fig. p.2| - Fig. |C.4| , observational data of ^He abundance, Y, obtained from metal- 
poor extragalactic Hnregion are seen as a function of metalicity (0/H, N/H, C/H) 
Fitting these data linearly, one can get the following fitting formulae; 



Y 
Y 



0.229 ± 0.004 + (1.3 ± 0.3) x 102(O/H), 
0.231 ± 0.003 + (2.8 ± 0.7) x 10^(N/H), 



(C.24) 
(C.25) 
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Figure C.2: He mass fraction Y vs. observed oxygen abundance. 
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Figure C.3: ^He mass fraction Y vs. observed nitrogen abundance. 
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Figure C.4: *^He mass fraction Y vs. observed carbon abundance. 
Y = 0.230 ±0.007+ (3.1 ± 2.2) X 102(C/H). 



(C.26) 



Prom these data, primordial abundance of ^He should be read off. Since the primordial 
^He in these metal-poor region is still contaminated by the one originating to stars, we take 
a limit that metalicity goes to zero. Then, the primordial abundance is estimated to be 






yp 


N 




C 





Yp = 0.229 d= 0.004, 
Yp = 0.231 ± 0.003, 



(C.27) 
(C.28) 
(C.29) 

Taking the observational uncertainty into account, we adopt the primordial ''He abundance 
as r„ = 0.23 ± 0.01. 



Next we will see the primordial abundance of ^Li. Since ^Li can be produced after BBN 
by cosmic ray spallation or some stellar processes (like novae outbursts), and it is easily 
destroyed at T > 2 x lO^K, we can not regard the present abundance of ^Li as primordial 
one. This means that we have to observe the old environment so as to see the primordial 
^Li abundance. For this purpose, observational data of ^Li abundance in metal-poor (and 
Tsuri ~ 5500K, as will be explained below) population II stars are usually used since they are 
expected to reflect the uncontaminated ^Li abundance. 
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Figure C.5: ^Li abundance [^Li] = 12 + log]^Q(''Li/H) in the metal-poor population II stars 
as a function of surface temperature Tgurf • The filled marks represent upper limits to the ''Li 
abundance. 



In Fig. |C.5| , we plot the ^Li abundance [''Li] = 12 + log^Q(^Li/H) in the most metal-poor 
population II stars as a function of surface temperature Tgurf- As one can see, "^Li-plateau" 
appears at Tsurf ~5500K. Using the data with Tsm-f > 5500K, the "plateau" value of [^Li] is 
estimated to be 



^Li 



plateau 



2.08 ±0.07. 



(C.30) 



We identify this value as the primordial abundance of ^Li. Notice that the data taken from 
stars with of Tgurf ^ 5500K are not appropriate for our purpose, because cool stars have thick 
surface convective zones which carry ^Li to deep region hot enough to burn it. Thus, we 
expect that ''Li abundances at the low surface temperature stars are not primordial one. 
Using eq. (|C.30|) , upperbound on the primordial ^Li abundance is given by 



^Li 



< 2.15. 



(C.31) 



Comparison with the theoretical predictions 

Now we are at the position to compare the observational data of light nuclei with the 
theoretical predictions. At first, we will constrain baryon-to-photon ratio tib from the data 
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of D, (D+^He) ^He, and ^Li. 

Let us begin with D and ^He. Lowerbound on the primordial abundance of D is given in 
eq.( C.21 ) and using that, r^B is constrained to be 



D : r^B < 6.8 X 10"^°. (C.32) 

On the other hand, upperbound on the primordial (D+^He) abundance, eq.( |C.23| ), requires 

D + ^He : r/B > 2.8 X 10"^°. (C.33) 

Theoretical prediction on the primordial abundance of ^He is considerably affected by 
an uncertainty of the neutron life time (r„ = 889.1 ± 2.1sec). In order to get a conservative 
constraint on the baryon-to-photon ratio rjB, we use the value = 891.2sec in deriving the 
lowerbound on rjB, and r„ = 887.0sec for upperbound. As a result, we obtain constraints on 
r]B from Yp = 0.23 ± 0.01 as 

^He : 1.5 x 10"^° <Vb< 3.3 x 10^^°. (C.34) 



Upperbound on the primordial ^Li abundance is given in eq.( |C3lD ; [^Li]p < 2.15. If we 



naively use this value, 7]b is constrained to be 1.9 x 10~^° ^ Vb ^ 3.3 x 10~^°. But from 
the uncertainties in the nuclear reaction rate, the predicted ^Li abundance is expected to 
be uncertain by ~40%. Assuming 40% residual uncertainty in the primordial abundance, 
bound on ?7b is found to be 

^Li : 1.6 X 10'^° <Vb< 4.0 x 10'^°. (C.35) 

Combining eq. (|C.32|) - eq. ( |C.35| ), allowed range of baryon-to-photon ratio rjs is given by 



D, ^He, ^He, ^Li : 2.8 x 10"^° < Vb < 3.3 x 10"^°, (C.36) 

i.e. riB ^ 3 X 10^^^ is predicted from BBN. 

Finally we will comment on the baryonic density of the present universe. By using 
eq. (|C.12|) , constraints on rjB (|C.36| ) becomes 

1.0 X 10-2 < QBh'^ < 1-2 X 10-2. (C.37) 

Therefore, baryonic dark matter {Qb ~ 1) conflicts with the theoretical predictions of BBN 
provided 0.5 ^ /i ^ 1. 



122 



Bibliography 



[1] J. Wess and B. Zumino, NucL Phys. B70 (1974) 39. 

[2] L. Maiani, Gif-sur-Yvette Summer School on Particle Physics, (Natl. Inst. Nucl. Phys. 
Part. Phys., Paris, 1979). 

[3] M.J.G. Veltman, Acta Phys. Pol. B12 (1981) 437. 

[4] P. Langacker and M. Luo, Phys. Rev. D44 (1991) 817. 

[5] U. Amaldi, W. de Boer and H. Fiirstenau, Phys. Lett. B260 (1991) 447. 

[6] S. Dimopoulos and H. Georgi, Nucl. Phys. B193 (1981) 150. 

[7] N. Sakai, Z. Phys. Cll (1981) 153. 

[8] E. Cremmer, S. Ferrara, L. Grardello and A. van Proeyen, Nucl. Phys. B212 (1983) 
413. 

[9] J. Elhs, C. Kounnas and D.V. Nanopoulos, Phys. Lett. B143 (1984) 410. 

[10] J. EUis, K. Enqvist and D.V. Nanopoulos, Phys. Lett. B147 (1984) 99. 

[11] S. Weinberg, Phys. Rev. Lett. 48 (1982) 1303. 

[12] L.M. Krauss, Nucl. Phys. B277 (1983) 556. 

[13] H. Pagels and J.R. Primack, Phys. Rev. Lett. 48 (1982) 223. 

[14] A.H. Guth, Phys. Rev. D23 (1981) 347. 

[15] J. Ellis, A.D. Linde and D.V. Nanopoulos, Phys. Lett. B118 (1982) 59. 

[16] T. Yanagida, Talk given at Physical Society of Japan Symposium on JLC, (Yamagata, 
October 1994). 



[17] H. Murayama, preprint LBL-36175, (LBL, 1994, |hep-ph/9410285D . 



[18] M.T. Grisaru, W. Siegel and M. Rocek, Nucl. Phys. B159 (1979) 429. 
[19] B.A. Campbell, S. Davidson, J. Ellis and K.A. Olive, Phys. Lett. B256 (1991) 457. 

123 



[20] H. Georgi and S.L. Glashow, Phys. Rev. Lett. 25 (1974) 438. 

[21] S. Kelly, J. L. Lopez and D. V. Nanopoulos, Phys. Lett. B247 (1992) 387. 

[22] V. Barger, M.S. Berger and P. Ohmann, Phys. Rev. D47 (1993) 1093. 

[23] P. Langacker and N. Polonsky, Phys. Rev. D49 (1994) 1454. 

[24] Particle Data Group, Phys. Rev. D50 (1994) 1173. 

[25] J. Gasser and H. Leutwyler, Phys. Rep. C87 (1982) 77. 

[26] J. Wess and J. Bagger, Supersymmetry and Supergragity, (Princeton University Press, 
1992). 

[27] J. Polonyi, preprint KFK-1977-93 (Budapest, 1977). 

[28] T. Banks, D.B. Kaplan and A.E. Nelson, Phys. Rev. D49 (1994) 779. 

[29] I. Joiclii, Y. Kawamura and M. Yamaguchi, preprint TU-462, (Tohoku University, 
1994, |hep-ph/ 94073851) . 

[30] G.D. Coughlan, W. Fischler, E.W. Kolb, S. Raby and G.G. Ross, Phys. Lett. B131 
(1983) 59. 

[31] P.R. Auvil and J.J. Brehm, Phys. Rev. 145 (1966) 1152. 
[32] P. Fayet, Phys. Lett. B175 (1986) 471. 

[33] D.A. Dicus, S. Nandi and J. Woodside, Phys. Lett. B258 (1991) 231. 

[34] D.V. Nanopoulos, K.A. Olive and M. Srednicki, Phys. Lett. B127 (1983) 30. 

[35] M.Yu. Khlopov and A.D. Linde, Phys. Lett. B138 (1984) 265. 

[36] J. Ellis, E. Kim and D.V. Nanopoulos, Phys. Lett. B145 (1984) 181. 

[37] R. Juszkiewicz, J. Silk and A. Stebbins, Phys. Lett. B158 (1985) 463. 

[38] J. Ellis, D.V. Nanopoulos and S. Sarkar, NucL Phys. B259 (1985) 175. 

[39] M. Kawasaki and K. Sato, Phys. Lett. B189 (1987) 23. 

[40] J. Ellis, G.B. Gelmini, J.L. Lopez, D.V. Nanopoulos and S. Sarker, Nucl. Phys. B373 
(1992) 399. 



[41] M. Kawasaki and T. Moroi, preprint TU-457, (Tohoku University, 1994, [hep- 
ph/940336^ ). 



124 



[42] T. Moroi, H. Murayama and M. Yamaguchi, Phys. Lett. B303 (1993) 289. 

[43] T. Moroi, M. Yamaguchi and T. Yanagida, Phys. Lett. B342 (1995) 105. 

[44] I. Affleck and M. Dine, Nud. Phys. B249 (1985) 361. 

[45] H. Murayama, H. Suzuki, T. Yanagida and J. Yokoyama, Phys. Rev. D50 (1994) 2356. 

[46] T. Moroi and T. Yanagida, Prog. Theor. Phys. 91 (1994) 1277. 



[47] M. Kawasaki and T. Moroi, preprint TU-463, (Tohoku University, 1994, [hep- 
ph/9408"32l| ). 



[48] W. Fischler, Phys. Lett. B332 (1994) 277. 

[49] T.P. Walker, G. Steigman, D.N. Schramm, K.A. Ohve and H.-S. Kang, Ap. J. 376 
(1991) 51. 

[50] R.D. Evans, The Atomic Nucleus (McGraw-Hill, New York, 1955). 

[51] R. Pfiffer, Z. Phys. 208 (1968) 129. 

[52] D.D. Paul, B.L. Herman, P. Mayer and D.L. Olson, Phys. Rev. Lett. 44 (1980) 44. 

[53] A.N. Gorbunov and A.T. Varfolomeev, Phys. Lett. 11 (1964) 137. 

[54] Yu.M. Arkatov, P.I. Vatset, V.I. Voloshchuk, V.A. Zolenko, I.M. Prokhorets and V.I. 
Chimil', Sov. J. Nucl. Phys. 19 (1974) 589. 

[55] J.D. Irish, R.G. Johnson, B.L. Berman, B.J. Thomas, K.G. McNeill and J.W. Jury, 
Can. J. Phys. 53 (1975) 802. 

[56] C.K. Malcolm, D.B. Webb, Y.M. Shin and D.M. Skopik, Phys. Lett. B47 (1973) 433. 

[57] L. Kawano, preprint FERMILAB-Pub-92/04-A (Fermilab., 1992). 

[58] J. Gratsias, R.J. Scherrer and D.N. Spergel, Phys. Lett. B262 (1991) 198. 

[59] M.H. Reno and D. Seckel, Phys. Rev. D37 (1988) 3441. 

[60] S. Dimopoulos, R. Esmailzadeh, L.J. Hall and G.D. Starkman, Ap. J. 330 (1988) 545. 

[61] S. Dimopoulos, R. Esmailzadeh, L.J. Hall and G.D. Starkman, NucL Phys. B311 
(1989) 699. 

[62] K. Hidaka, Phys. Rev. D44 (1991) 927. 

[63] J. McDonald, K.A. Olive and M. Srednicki, Phys. Lett. B283 (1992) 80. 



125 



[64] S. Mizuta and M. Yamaguchi, Phys. Lett. B298 (1993) 120. 

[65] M. Drees and M.M. Nojiri, Phys. Rev. D47 (1993) 376. 

[66] J. Ellis, C. Kounnas and D.V. Nanopoulos, NucL Phys. B247 (1984) 373. 

[67] K. Inoue, M. Kawasaki, M. Yamaguchi and T. Yanagida, Phys. Rev. D45 (1992) 328. 

[68] S. Kelley, J.L. Lopez, D.V. Nanopoulos, H. Pois and K. Yuan, Nucl. Phys. B398 (1993) 
3. 

[69] E. Witten, Phys. Lett. B155 (1985) 151. 

[70] M. Matsumoto, J. Arafune, H. Tanaka and K. Shiraishi, Phys. Rev. D46 (1992) 3966. 

[71] J.L. Lopez, D.V. Nanopoulos and A. Zichichi, Phys. Lett. B291 (1992) 255. 

[72] J.L. Lopez, D.V. Nanopoulos and H. Pois, Phys. Rev. D47 (1993) 2468. 

[73] P.J. Kernan and L.M. Krauss, Phys. Rev. Lett. 72 (1994) 3309. 

[74] M. Fukugita and N. Sakai, Phys. Lett. B114 (1982) 23. 

[75] A. A. Zdziarski, Ap. J. 335 (1988) 786. 

[76] R. Svensson and A. A. Zdziarski, Ap. J. 349 (1990) 415. 

[77] A.D. Linde, Phys. Lett. B129 (1983) 177. 

[78] D.S. Salopek, Phys. Rev. Lett. 69 (1992) 3602. 

[79] E.W. Kolb and M.S. Turner, The Early Universe (Addison- Wesley, 1990). 

[80] B.W. Lee and S. Weinberg, Phys. Rev. Lett. 39 (1977) 165. 

[81] N. Terasawa, M. Kawasaki and K. Sato, Nucl. Phys. B302 (1988) 697. 

[82] J. Preskill, M.B. Wise and F. Wilczek, Phys. Lett. B120 (1983) 127. 

[83] L.F. Abbott and P. Sikivie, Phys. Lett. B120 (1983) 133. 

[84] M. Dine and W. Fischler, Phys. Lett. B120 (1983) 137. 

[85] F.A. Agaronyan, A.M. Atyan and A.N. Najapetyan, Astrofizika 19 (1983) 323. 
[86] L.C. Maximon, J. Res. NBS 72(B) (1968) 79. 

[87] V.B. Berestetskii, E.M. Lifshitz and L.P. Pitaevskii, Quantum Electrodynamics (Perg- 
amon press, Oxford, 1982). 



126 



[88] F.C. Jones, Phys. Rev. 167 (1968) 1159. 



[89] M. Kawasaki and T. Moroi, preprint TU-474, (Tohoku University, 1994, |astro- 
ph/941 20551 ). 



127 



